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ABSTRACT 

We have mapped the molecular gas content in the host galaxy of the strongly lensed high redshift 
quasar APM 08279+5255 (z = 3.911) with the Very Large Array at 0.3" resolution. The CO(J=1--^0) 
emission is clearly resolved in our maps. The C0( J=1^0) line luminosity derived from these maps is 
in good agreement with a previous single-dish measurement. In contrast to previous interferometer- 
based studies, we find that the full molecular gas reservoir is situated in two compact peaks separated 
by ;<0.4". Our observations reveal, for the first time, that the emission from cold molecular gas 
is virtually cospatial with the optical/near-infrared continuum emission of the central AGN in this 
source. This striking similarity in morphology indicates that the molecular gas is situated in a compact 
region close to the AGN. Based on the high resolution CO maps, we present a revised model for the 
gravitational lensing in this system, which indicates that the molecular gas emission is magnified by 
only a factor of 4 (in contrast to previously suggested factors of 100). This model suggests that the 
CO is situated in a circumnuclear disk of '^550 pc radius that is possibly seen at an inclination of 
<25°, i.e., relatively close to face-on. From the CO luminosity, we derive a molecular gas mass of 
Afgas=1.3x 10^^ Mq for this galaxy. From the CO structure and linewidth, we derive a dynamical 
mass of Mdyn sin^z=4.0x 10^° M0. Based on a revised mass estimate for the central black hole of 
il/BH=2.3xlO^° Mq and the results of our molecular line study, we find that the mass of the stellar 
bulge of APM 08279+5255 falls short of the local MBH-^buige relationship of nearby galaxies by more 
than an order of magnitude, lending support to recent suggestions that this relation may evolve with 
cosmic time and/or change toward the high mass end. 

Subject headings: galaxies: active — galaxies: starburst — galaxies: formation — galaxies: high- 
redshift — cosmology: observations — radio lines: galaxies 



1. INTRODUCTION 

Several different populations of distant galaxies have 
been detected to date, out to a (spectroscopically con- 
firmed) redshift of z=7.0 (see Ellis 2007 for a review). 
It is of fundamental importance to study such young 
galaxies in great detail in order to constrain their nu- 
clear, stellar, and gaseous constituents and their phys- 
ical properties and chemical abundances. Understand- 
ing these different characteristics of high redshift galaxy 
populations and their progression through cosmic times 
is vital to develop a unified picture of galaxy formation 
and evolution. One key aspect in these high-z galaxy 
studies are sensitive, high-resolution radio observations 
of the molecular gas phase, i.e., the raw material that 
fuels star formation. 

Since its initial discovery at z> 2 more than a decade 
ago (Brown & Vanden BoutEMl Solomon et al. ll992f ). 
observations of spectral line emission from interstellar 
molecular gas in distant galaxies have revolutionized our 
understanding of some of the most luminous objects that 
populate the early universe. This is due to the fact that 
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the physical state of the molecular interstellar medium 
(ISM) plays a critical role in the evolution of a galaxy. 
The total amount of molecular gas in a galaxy determines 
for how long starburst activity can be maintained, while 
its temperature and density are a direct measure for the 
conditions under which star f ormation can occur (see 
Solomon & Vanden Bout l2005l for a review). Also, there 
is growing evidence that the observed relationships be- 
tween black hole mass and galaxy bulg e veloc ity disper- 
sion (MBH-o-buigo; Ferrarese & Merritt [2000t Gebhardt 
et al. I200Q) . black hole mass and host bulge concentra- 



tion (i .e., Se rsic index; Graham et al. 120011 : Graham & 



Driver [2b07f) . and, ultimately, black hole m ass and bulge 
mass (MBH-Afbuige; Magorrian et al. ll998l : Haring & Rix 
may be a consequence of an active galactic nucleus 
(AGN) feedback mechanism acting on its surrounding 
mate rial, a t least in the mo st lum inous systems (Silk & 
Rees ll998t Di Matteo et al. l2005| ). Such an AGN-driven 
wind will also interact with the molecular gas, the ma- 
terial which will eventually form the stellar bulges and 
disks in galaxies at high redshift, but also feed the active 
nucleus itself. Observations of the dynamical structure 
and distribution of molecular gas in distant galaxies may 
even reveal the initial cause of event (e.g., mergers) for 
both star formation and AGN activity, an important test 
for recent cosmological models (e.g., Springel et al. 200l; 
Narayanan et al. 2006). 

The flux magnification provided by gravitational lens- 
ing has been facilitated in studies of high redshift molec- 
ular gas emission since early on, leading to the detection 
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of a num ber of quasar host g alaxies (e.g., Barvainis et 
al. Il994t Downes et al. Il999l ) and submillimeter galax- 
ies (SMGs; e.g., Frayer et al. Il998t Greve et al. I2005D . 
but also substantially fainter popula tions s uch as Lyman- 
break galaxies (LBGs; Baker et al. [2004'; Coppin et al. 
l2007t ). Also, some sources have large enough lens im- 
age separations to show structu re eve n at only m oderat e 
angular resolution (Sheth et al. 12004 Kneib et al. l2005D . 

High angular resolution observations of molecular gas 
in high-z galaxies are, to date, exclusively obtained 
in the rotational transitions of carbon monoxide (CO). 
In the p ast few years, a number of SMGs (Ge nzel et 
al. 120031 : Downes & Solo mon jl OOl Tacconi et al. 120061 ) 
and quasars (AUoin et al. 1997; at z>2 have been stud- 
ied at a linear resolution of up to 5kpc at the target 
redshifts (lkpc=0.12" at z=2). In some cases, gravita- 
tional lensing aids in zooming in further on these sys- 
tems. However, the only telescope that currently al- 
lows to attain resolutions of sub-kpc to Ikpc scale at 
even higher redshifts {z>4, where lkpc>0.14") is the 
NRAO's Very Large Array (VLA)^. First sub-arcsecond 
resolution CO imaging with the VLA has been obtained 
toward th e dist ant quasars BR 1202-0725 (z=4.69; Car- 
ilh et_aL UMl), APM 08279-1-5255 (z=3.91; Le wis et 
aL l2002aD . PS S J2322 -H944 (z=4.12; CariUi et al. [20031 : 
Riechers eiaL[20Mi), SDSS J1148-I-5251 (z=6.42; Wal- 
ter et al. j2004f ). and BRI 1335-0417 (z=4.41; Riechers et 



TABLE 1 

New (2005) and archival (1998-2001) 
vla observations summary. 



al. I2008bl) 



We here report on new, more sensitive high-resolution 
observations of CO(J=1^0) emission towards the 
strongly lensed quasar APM 08279+5255 (z = 3.911). 
APM 08279-1-5255 was already^ imaged during the Palo- 
mar Sky Survey (PSS) in 1953, a decade before the 
discovery of th e first quasar was reported in the liter- 
ature f Schmidt [1961) . About 45 years later, it was 'of- 
ficially' discovered in an automatic plate measuring fa- 
cility (APM) survey for distant cool carbon stars, and 
identified as a 15"^ magnitude, radio-quiet broad ab- 
sorption line (BAL) quasar with an IRAS Faint Source 
Catalog (FSC) counterpart at a redshift of Zopt=3.87 
(Irwin et al. 1998). APM 08279-1-5255 was found to 
have an unprecedented apparent bolometric luminosity 
of Lboi=7x 10^^ Lq, about 20% of wh ich is emitted in 
the (far-) infrared (Lewis et al. Il998[ ). This extreme 
value was found early on to be due to strong gravita- 
tional lensing (Irwin et al. 1998), producing three close 
(maximum separation <0.4"), almost co llinea r images in 
th e opti cal/near-infrared (Ibata et al. 119991 : Egami et 
al. |2000| ). Due to the lack of a detection of the lensing 
galaxy, and an indication for significant microlensing ef- 
fects, the true nature of the gravitational lens configura- 
tion (and thus the magnification factor hl) in this system 
remai ns subje ct to debate (e.g., Egami et al. 120001 : Lewis 
et al. I2002£j ) despite the fact that APM 08279-^5255 
is one of the best-studied sources in the distant uni- 
verse. In addition to its enormous continuum brightness 
at basically every astronomically relevant wavelength, 
APM 08279-^5255 is also one of the brightest CO sources 

^ The Very Large Array is a facility of the National Radio As- 
tronomy Observatory, operated by Associated Universities, Inc., 
under a cooperative agreement with the National Science Founda- 
tion. 

Due to the source's enormous optical brightness, much older 
images may exist. 



Date 


Config. 


Duration 


Band 


Type 




Ref. 


2005 Apr 07 


B 


6.75 h 


K 


CO(J= 


1- 


-0) 


1 


2005 Apr 08 


B 


6.0h 


K 


CO(J= 


1- 


-0) 


1 


2005 Apr 09 


B 


5.25 h 


K 


CO(J= 


1- 


-0) 


1 


2005 Apr 10 


B 


6.0h 


K 


CO(J= 


1- 


-0) 


1 


2005 Apr 18 


B 


7.5 h 


K 


CO J= 


1- 


-0) 


1 


2005 May 01 


B 


4. Oh 


C 


Continuum 


1 


2005 Sep 12 


C 


8.25 h 


K 


CO(J= 


1- 


-0) 


1 


2005 Sep 20 


C 


7.5 h 


K 


CO(J= 


1- 


-0) 


1 


2005 Dec 02 


D 


4.0h 


K 


CO(J= 


1- 


-0) 


1 


1998 Jun 18 


BnA 


3.0h 


X 


Continuum 


2 


1998 Jul 28 


B 


2.5 h 


X 


Continuum 


1 






2.0h 


u 


Continuum 


1 


1998 Aug 11 


B 


3.0h 


X 


Continuum 


1 






3.5 h 


u 


Continuum 


1 


2000 Apr 24 


C 


13. Oh 


Q 


CO(J= 


2- 


-1) 


3 


2000 Sep 16 


D 


4.0h 


Q 


CO J= 


2- 


-1) 


1 


2001 Jan 21 


A 


2.25 h 


L 


Continuum 


4 


2001 Oct 20 


D 


9.0h 


Q 


CO(J= 


2- 


-1) 


1 



Re feren ces. — [1] this work; [2] Ibata ct al.j2999,; [3] Papadopoulos 
ct al. 120011 : [4] Ivison [2006l . 

Note. — Maximum baselines of VLA configurations are: A -36.4 km, 
B - 11.4km, C - 3.4km, D - 1.03km. BnA is a hybrid configuration 
where the antennas on the east and west arms of the array are in the 
shorter of the two denoted configurations, and the north arm is in the 
more extended configuration. 



at high redshift (Downes et al. 11999 ]^ It also exhibits 

unusually bright HCN(Wagg et al. l2005f) and HCO+ 
(Garcia- Burillo ct al. l2006l ) emiss ion, a nd was recently 
also detected in C i (Wagg et al. l2006f ) emission. The 
properties of its molecular gas content have recently been 
modeled based on observations of the C0(J=1— >0) up 
to CO(J=11^10) rotational ladder ('CO spectral line 
energy distribution', or SLED) with the NRAO Green 
Bank Telescope (GBT)^ and the IRAMSOm telescope 
(Riechers et al. I2OO6I : Weifi et al. 1.20071). We use a 
standard concordance cosmology throughout, with Hq — 
71 k ms-\Mpc- \ Om = 0.27, and Qa = 0.73 (Spergel 
etal. [200a [2001 . 

2. OBSERVATIONS 

We here report new, sensitive high-resolution ob- 
servations of the CO(J=1^0) emission line in 
APM 08279-1-5255. To better constrain the molecular gas 
properties in this source, we also re-derive the strength 
of the C0(J=2— >1) emission line, based on more sen- 
sitive (unpublished archival) obs ervati ons than reported 
previously (Papadopoulos et al. l2001f ). In this context, 
we also present a search for the CN(A^=1— >0) emission 
line in this galaxy, a tracer of the dense fraction of the 
molecular gas phase (see Riechers et al. [2006 for more de- 
tails on this data). To further investigate the AGN and 
starburst environments in this source, we augment our 
analysis of the molecular gas properties with a study of 
the radio continuum properties, based in part on archival 
observations. A fraction of the archival continuum data 
was u sed in previ ous studies of this source (Ibata et al. 
119991: Ivison|2006t see Tab.[T]for details). 

* The Green Bank Telescope is a facility of the National Radio 
Astronomy Observatory, operated by Associated Universities, Inc., 
under a cooperative agreement with the National Science Founda- 
tion. 



High-Resolution CO Imaging of APM 08279+5255 
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2.1. VLA Data 
2.1.1. CO (J =1^0) and 23.5 GHz continuum [new] 

We observed the CO(J=1^0) transition hne (^'rcst = 
115.2712018 GHz) towards APM 08279+5255 using the 
VLA in B, C, and D configuration between 2005 April 
07 and 2005 December 02. At the target redshift of 3.911, 
the line is redshifted to 23.472GHz (12.77mm, K band). 
The total on-sky integration time amounts to 51.25 hr 
(see Tab. [1] for details). Observations were perfor med in 
fast-switching mode (see, e.g., Carilli & Holdawav 1 1999( 1 
using the nearby source 08248+55527 for secondary am- 
plitude and phase calibration. Observations were carried 
out under excellent weather conditions with 25 anten- 
nas. The phase stability in all runs was excellent (typi- 
cally <15° for the longest baselines). The phase coher- 
ence was determined by imaging the calibrator source 
08087+49506 with the same calibration cycle as that 
used for the target source, and found to be high dur- 
ing all runs. For primary flux calibration, 3C286 was 
observed during each run. 

Two 25 MHz intermediate frequencies (IFs) with seven 
3.125 MHz channels each were observed in 'spectral line' 
mode centered at the C0(J=1— >0) line frequency, lead- 
ing to an effective bandwidth of 43.75 MHz (correspond- 
ing to 558 kms"^ at 23.5 GHz; 'ON' channels). Two 
50 MHz (corresponding to 638kms~i at 23.5GHz) IFs 
were observed in continuum mode at 23.3649 GHz and 
23.5649 GHz (±100 MHz offset from the line frequency) 
to measure the source's continuum at 12.8 mm ('OFF' 
channels). The continuum was observed for one third of 
the total time to attain the same rms as in the combined 
line channels. 

For data reduction and analysis, the AIVS'^ package 
was used. The two continuum channels were concate- 
nated in the uv/visibility plane. The CO(J=1^0) line 
image was generated by subtracting a CLEAN compo- 
nent model of the continuum emission from the visibility 
data. All data were mapped using the CLEAN algo- 
rithm. The velocity-integrated maps shown in Fig. [1] 
are imaged using robust 1 weighting and combining B 
and C array data, leading to a synthesized clean beam 
size of 0.30 "x 0.29" in the CO(J=1^0) line map and 
0.32 "x 0.29" in the continuum map (note that D array 
data were omitted in these maps to improve the size of 
the synthesized beam) . The final rms over the full band- 
width of 43.75MHz (558 kms'^) in the CO(J=1^0) 
map is 16/iJy beam"^, and the rms over the 100 MHz 
continuum bandpass is 13 ^Jy beam~^. In Fig. [21 four 
velocity channel maps (9.375 MHz, or 120 kms~^ each) 
of the central 37.5 MHz (480 kms^^) of the CO(J=1^0) 
line are shown. To maximize the sensitivity in individual 
channels the data were imaged using natural weighting 
and combining B and C array data, achieving a resolu- 
tion of 0.34"x0.33". Banning smoothing was applied 
to the image cube, resulting in an rms of 27/iJybeam^^ 
per channel. The velocity-integrated C0(J=1— >0) and 
23.5 GHz continuum maps of the combined B, C and 
D array data shown in Fig. [3] and Fig. |4l top right, 
are imaged using natural weighting and a Gaussian ta- 
per, leading to synthesized beams of 1.19 "x 1.14" and 
1.19"xl.l0", and an rms of 19^Jy beam^^ for both 

S www.aoc.nrao.edu/aips/ 



TABLE 2 

CO/CN Line Peak and Continuum Fluxes. 



Transition /Frequency 


[^Jy] 


Configuration 


CO(J=1^0) 


285+27'' 


BCD 


CO(J=2^1) 


1637+65'' 


CD 


CN(Af=1^0) 


<195 


[GBT] 


1.4 GHz 

4.5 GHz 
8.4 GHz 
14.9 GHz 
23.5 GHz 

46.9 GHz (estimated)^ 


1160+33 
551+40 
446+20 
303+93^^ 
376+19 
405+310 


A 
B 

BnA/B 
B 
BCD 



Note. — Fluxes are derived from convolved maps for 
all datasets with linear resolutions higher than 1.2", and 
agree with those derived from Gaussian fitting to the 
unconvolved data within 10%. 

^ Peak flux from Gaussian fitting to the unsmoothed 14 
channel spectrum. 

^ Not corrected for the estimated 46.9 GHz continuum 

contribution. 

^ Tentative detection. 

^ Estimated based on the FIR-to-radio SED. The lower 
limit is given by the model of the FIR-only part of the 
continuum, and the upper limit is estimated by linear 
interpolation (in the log— log plane) of the 23 GHz and 
94 GHz continuum measurements. 

maps. 

2.1.2. CO(J=2^1) [arctiival] 

Spectral line observations of the C0(J=2--i-l) 
transition line (i^rest — 230.53799 GHz) towards 
APM 08279+5255 were carried out between 2000 Apr 24 
and 2001 Oct 20 using the VLA in C and D configura- 
tion. The C arr ay data was published by Papadopou- 
los et al. ()2001f) . We here re-analyze this data, and 
add in shorter (unpublished) D array spacings to im- 
prove the sensitivity of the data. At the target red- 
shift of 3.911, the C0(J=2^1) line is redshifted to 
46.943GHz (6.39mm, Q band). The total on-sky in- 
tegration time amounts to 26 hr (see Tab. [T] for details). 
Observations were performed in fast-switching mode us- 
ing the nearby source 08248+55527 for secondary am- 
plitude and phase calibration. Observations were car- 
ried out under excellent weather conditions with up to 
27 antennas. For primary flux calibration, 3C 48 and 
3C 286 were observed. Observations were set up in quasi- 
continuum mode; two 50 MHz IFs were observed simul- 
taneously at 46.9149 GHz and 46.9649 GHz to cover the 
central '^640 km s""'^ of the spectral line. 

Data were mapped using the CLEAN algorithm at 
natural weighting; this results in a synthesized beam 
of l.ll"x0.93" and an rms of 65/iJy beam~^ (no map 
shown here). 

2.1.3. 1.4 GHz continuum [archival] 

Continuum observations at 1.400 GHz (L band) to- 
wards APM 08279+5255 were carried out on 2001 Jan 21 
using the VLA in A config uration. This data was recently 
published by Ivison (|2006t no map shown). The total on- 
sky integration time amounts to 2.25 hr (see Tab. [1] for 
details). Observations were performed in spectral line 
mode using the nearby source 0824+558 for secondary 
amplitude and phase calibration. Observations were set 
up using two IFs and seven 3.125 MHz channels per IF, 
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Fig. 1. — High-resolution VLA maps of CO(J=l— >0) + continuum (left), CO(J=l— »0) only (middle), and 23.5 GHz continuum (right) 
emission towards APM 08279+5255 at a resolution of 0.30/0.30/0.32 " xO.29 " (robust 1 weighting; left to right; beam size indicated in the 
bottom left corner of each panel). Contours are shown in steps of the rms noise, starting at ±2(t (Icr = 16/16/13 /ijy beam~^). Left: 
Image of the CO(J=1^0) and continuum emission. The source is clearly resolved, and detected at a peak strength of 15cr. This map is 
integrated over the central 558kms~^ (43.75 MHz) of the CO(J=1^0) emission. Middle: Image of the CO(J=l— >0) emission. Continuum 
emission was subtracted. The source is clearly resolved into two CO peaks (8cr and 5(t signal-to- noise), as indicated by the plus signs. This 
continuum-subtracted map is integrated over the same range as in left the panel. Right: Image of the 23.5 GHz continuum emission (peak: 
llcr). The emission appears resolved and peaks centrally between the CO(J=l— >0) peaks (as indicated by the plus signs). This continuum 
map is integrated over a bandwidth of 100 MHz. 



leading to an effective bandwidth of 43.75 MHz (corre- 
sponding to ^9400 km s""'^ at 1.4 GHz). Observations 
were carried out under good weather conditions with 27 
antennas. For primary flux cahbration, 3C 48 was ob- 
served. 

Data were mapped using the CLEAN algorithm at ro- 
bust weighting; this results in a synthesized beam of 
1.42"xl.l5" and an rms of 33/iJy beam"^ (see Fig. gl 
top left). 

2.1.4. 4-5 GHz continuum [new] 

Continuum observations at 4.5276 GHz (C band) to- 
wards APM 08279-^5255 were carried out on 2005 May 
01 using the VLA in B configuration. The total on-sky 
integration time amounts to 4hr (see Tab. [T] for details). 
Observations were performed in spectral line mode using 
the nearby source 08248-1-55527 for secondary amplitude 
and phase calibration. Observations at 4.5276 GHz were 
set up using thirty-one 390.625 kHz channels, leading to 
an effective bandwidth of 12. 109375 MHz (corresponding 
to 802kms^i at 4.5276 GHz). Observations were carried 
out under excellent weather conditions with 25 antennas. 
The phase stability in this run was excellent. For primary 
flux calibration, 3C 286 was observed during each run. 

Data were mapped using the CLEAN algorithm at 
natural weighting; this results in a synthesized beam of 
1.55 "x 1.42" and an rms of 40 ^Jy beam~^ (see Fig. |31 
bottom left). 

2.1.5. 8.4 GHz continuum [archival] 

Continuum observations at 8.4 GHz (X band) towards 
APM 08279+5255 were carried out between 1998 Jun 18 
and Aug 11 using the VLA in BnA and B configuration. 
The BnA array data was published by Ibata et al. (,1999 ). 
We here re-analyze this data, and add in shorter B array 
spacings to improve the sensitivity and image quality of 
the data. The total on-sky integration time amounts to 
8.5 hr (see Tab. [1] for details). Observations were per- 
formed in quasi-continuum mode (two 50 MHz IFs were 
observed simultaneously at 8.4351 GHz and 8.4851 GHz) 



using the nearby sources 0749+540 and 0820+560 for sec- 
ondary amplitude and phase calibration. Observations 
were carried out under good weather conditions with 27 
antennas. For primary flux calibration, 3C 48 and 3C 147 
were observed. 

Data were mapped using the CLEAN algorithm. Imag- 
ing the data with natural weighting results in a syn- 
thesized beam of 0.73 "x 0.42" and an rms of 10/iJy 
beam~^ (see Fig. SI bottom right). Imaging the data 
with uniform weighting results in a synthesized beam of 
0.63 "x 0.22" and an rms of 20/xJy beam"! (see Fig. [5]). 

2.1.6. 14-9 GHz continuum [archival] 

Continuum observations at 14.9 GHz (U band) towards 
APM 08279+5255 were carried out on 1998 Jul 28 and 
Aug 1 1 using the VLA in B configuration. The total on- 
sky integration time amounts to 5.5 hr (see Tab.[T]for de- 
tails). Observations were performed in quasi-continuum 
mode (two 50 MHz IFs were observed simultaneously at 
14.9149 GHz and 14.9649 GHz) using the nearby source 
0749+540 for secondary amplitude and phase calibration. 
Observations were carried out under good weather con- 
ditions with 27 antennas. For primary flux calibration, 
3C 48 was observed. 

Data were mapped using the CLEAN algorithm at 
natural weighting; this results in a synthesized beam 
of 0.51 "xO. 44" and an rms of 44/iJy beam^^ (no map 
shown here). 

2.2. GET Data 



CN(N=1^0) 

wide-ba nd h igh-z molecular 
et al. [2006h . we observed 
i^=5/2^3/2, F=l/2^1/2, 
hyper-fine structure transi- 
of the CN(iV=1^0, u=0, J=3/2^1/2) line 
= 113.4881420 - 113.5089340 GHz) towards 
APM 08279+5255 with the GBT between 2004 Decem- 
ber 15 and 19 for a total of 22.5 hr (15 hr on source). At 
the target redshift of 3.911, these lines are redshifted 
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Fig. 2. — Channel maps of the CO(J=l— >0) emission (same region is shown as in Fig. [T}. The beam size (0.34 " xO.33 "; natural 
weighting) is shown in the bottom left corner. The plus signs indicate the same positions as in Fig. [T] Contours are shown in steps of 
the rms noise, starting at ±2fT (Ic = 27/iJy beam~^). One channel width is 9.375 MHz, or 120kms~^ (frequencies increase with channel 
number and are shown at 23451.539, 23460.914, 23470.289, and 23479.664 MHz). 



to 23.109-23. 113GHz (12.97mm, K band), and can be 
observed simultaneously (in-ban d) wit h C0(J=1— >0) 
at the GET (see Riechers et al. l2006l for a more de- 
tailed description and results on the CO observations). 
Observations were performed in ON-OFF position 
switching mode using the nearby source J0753-I-538 
for secondary amplitude and pointing calibration, 
correcting for the atmospheric opacity. The pointing 
accuracy was typically 3". Observations were carried 
out under excellent weather conditions. For primary 
flux calibration, 3C 147 and 3C 286 were observed during 
each run. The dual-beam, dual-polarization 18-26 GHz 
receiver was used for all observations. The beam size at 
our observing frequency is 32" ('^230kpc at the source's 
redshift), i.e., much larger than our target. 

Two 800MHz IFs with 2048 390.625kHz (5 kms'^) 
channels and two orthogonal polarizations each were ob- 
served simultaneously. The first IF was centered on the 
redshifted CO(J=1^0) frequency at 23.4663 GHz, and 
the second IF was tuned 700 MHz bluewards. Therefore, 
the observed CN(iV=1^0) fine structure transition lines 
are located at the lower end of the first IF. 

For data reduction and analysis, the AIPS++^" and 
GET IDL^^ packages were used. Final combina- 
tion, binning, and baseline subtraction were done with 
GILDAS/CLASS^^. An rms of 65 ^Jy was attained in 
these observations. 

3. RESULTS 

3.1. CO(l-O) Line and 23.5 GHz Continuum Emission 

We have resolved the emission from the CO(J=1^0) 
transition line towards APM 08279-^5255 {z = 3.91), 
as well as the continuum emission at the same fre- 
quency (Fig. [T|). The left panel of Fig. [1] shows the 
C0( J=l— >0) and continuum emission, the middle panel 
the continuum-subtracted C0(J=1— >0) emission, and 
the right panel the 23.5 GHz continuum emission. The 
continuum-subtracted CO(J=1-^0) hue emission (Fig. 
[T] middle) is clearly resolved into two peaks (north-east 
and south-west). From 2-dimensional Gaussian fitting 
we find a line peak flux density of 131 ± 16/xJy beam"^ 
for the north-eastern CO peak at a=08h31'"41".708 ± 
0^001, J=-|-52°45'17".58 ± 0".02 and a line peak flux 
density of 87 ± 16/xJy beam^^ for the south-western CO 

http;//www. aips2.nrao.edu. 
http:/ /gbtidl. sourceforge.net. 
12 http://www.iram.fr/IRAMFR/GILDAS/. 



peak at a=08''31"41".688 ± 0".004, (5=-f 52°45'17".34 ± 
0".05 (astrometric errors are derived from the flt; Ta- 
bles [2] and [3|) . This corresponds to a peak flux ratio of 
6(SW,NE)==0.66 ± 0.15. The separation between both 
peaks is d(SW,NE)=0.31" ± 0.05". 

The continuum emission (Fig. [U right) is also 
clearly resolved and has about the same extent as the 
CO(J=1-^0) line emission, but is not resolved into two 
individual peaks. The emission rather peaks centrally 
between the two CO peaks indicated by the plus signs. 
For the continuum emission at this resolution, we derive 
a peak flux density of 147 ± 12 /^Jy beam~^ by fltting a 
single two-dimensional Gaussian. 

In Fig. [21 four (continuum-subtracted) channel maps 
(9.375MHz, or 120 kms"! each) of the central 37.5MHz 
(480 kms-i) of the CO(J=1^0) line are shown. The 
line is resolved into the two CO(J=1^0) peaks defined 
above in the central two channels, and the decline of the 
line intensity towards the line wings is clearly visible in 
the outer channels, as expected. Due to the slight asym- 
metry of the CO(J=1^0) line with respect to our central 
tuning frequency, only little emission is detected in the 
blueward outer channel, the lack of extended emission in 
this channel is therefore ascribed to low signal-to-noise 
rather than different structure relative to the other chan- 
nels. However, there is tentative evidence for a spatially 
resolved velocity gradient between the velocity channels 
along an axis close to the major lensing axis. From the 
central line channels, we derive a line peak flux density 
of 185 ± 27/iJy beam~^ for the north-eastern CO peak 
and a line peak flux density of 101 ± 27/xJy beam"^ for 
the south-western CO peak. 

To derive the fuU CO(J=1^0) and 23.5 GHz con- 
tinuum flux of APM 08279-1-5255, we have imaged the 
whole E, C, and D array dataset with natural weight- 
ing, but tapered to a linear resolution of ^-^1.2" (Fig. [3] 
and Fig. IH top right). At this resolution, the source 
appears at best marginally resolved, and we derive that 
we outresolve at most 10% of the integrated flux. From 
the integrated maps, we derive an CO(J=1^0) peak 
flux density of 210 ± 19 /iJy beam^^, and a contin- 
uum peak flux density of 372 ± 19/iJy beam^^. The 
complete beam-corrected source fluxes (derived from 2- 
dimensional Gaussian fitting to the source) in both maps 
agree with the peak fluxes within the errors. It also 
agrees well with the beam-corrected fluxes in the full res- 
olution dataset within the error bars. From the central 
line channels, we derive a CO(J=1-^0) line peak flux 
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Fig. 3. — Map of the CO(J=1-+0) emission as shown in Fig. [l] 
(middle), but including D array short spacings, and tapered to a 
resolution of 1.19"xl.l4" (natural weighting; beam size indicated 
in the bottom left corner). Contours are shown in steps of the 
rms noise, starting at ±3cr (la = 19 /ijy beam~^; llcr peak). The 
plus signs indicate the same positions as in Fig. [T] The emission 
appears practically unresolved, and there is no indication of any 
extended emission. 

density of 285 ± 27/xJy beam"-'^. 

This peak flux agrees well with that derived from 
single-dish observations. Figure [6] shows the central part 
of the C0( J=l— >0) spectrum we have obtained with the 
GBT (Riechers et al. |2006). This spectrum has a res- 
olution of 75 kms~^ at an rms of 65/xJy per channel. 
The solid black line is a Gaussian fit to that data. Over- 
plotted on this spectrum, we show the flux of the four 
120 kms~^ channels of Fig. [2] (but for the tapered BCD 
array dataset) . The error bars indicate the rms per chan- 
nel (40/iJy). This clearly demonstrates that, within the 
error bars, our interferometric observations recover the 
same amount of CO emission as the single-dish observa- 
tions obtained with a ^32" observing beam. 

3.2. C0(2~l) Line Emission 

The C0(J=2-^1) transition line towards 
APM 08279-^5255 {z = 3.91) is also det ected (see 
also discovery paper by Papadopoulos et al. l2001h . To 
extract the full C0(J=2-^1) line flux, we convolved 
the map with a circular beam of 1.2" diameter. From 
2-dimensional Gaussian fitting, we then derive a line 
peak flux density of 1637 ± 65/xJy beam^^. 

We note that during none of the (archival) 
C0(J=2^1) observations, the continuum at 46.9 G Hz 
was measured. However, Papadopoulos et al. ()200lD re- 
port a limit of 3cr=0.9mJy on the continuum emission 
at 43.3 GHz. This means that up to about half of the 
C0( J=2— >1) emission reported above may be attributed 
to 7mm continuum emission rather than C0(J=2— >1) 
emission (see Sect. 4.3 for further discussion). 

3.3. The Search for CN(l-O) Line Emission 

To measure the total CO luminosity, we have ob- 
tained wide-band CO(J=:1^0) observations toward 
APM 08279-1-5255 with the GBT (see Riechers et al. 



|2006| for details). The bandwidth of these observa- 
tions was high enough to also search for redshifted 
CN(A^=1— >0) emission within the same setup (see 
Fig. [7|). The strongest hyper-flne structure compo- 
nents of the CN(A^=1^0, w=0, J=3/2^1/2) multi- 
plet, i.e. the i^=3/2-+l/2, i^=5/2^3/2, F=l/2^1/2, 
and F=3/2— >3/2 transitions, have rest frequen- 
cies of 113.4881420, 113.4909850, 113.4996430, and 
113.5089340 GHz. At z = 3.911, this corresponds to 
23.109, 23.110, 23.111, and 23.113 GHz. Accounting for 
the uncertainty in redshift (Riechers et al. ,2006) , the 
dashed lines in Fig. [7] indicate the region where these 
transitions are expected to peak. We note that, if the 
CN emission comes from a similar region than the CO 
emission, these components are expected to be blended 
due to their large linewidth (>500 kms ~\ or 39 MHz 
FWHM, Riechers et al. l200l Weifi et al. 120071) . At an 
rms of 65/LtJy, we do not detect any emission. We thus 
set a 3(7 limit of 195 /iJy on the peak strength of the 
CN(Af=1^0) emission in APM 08279-1-5255. 

3.4. Continuum Emission at L4 GHz 

Using the VLA in A array, 1.4 GHz continuum emission 
from APM 08279-1-5255 was dete cted by Ivison (2006). 
From this analysis, Ivison l)2006f ) derived a continuum 
flux density of 3.05 ± 0.07 mJy. APM 08279-^5255 was 
previously detecte d at 1 .4 GHz in the VLA-FIRST sur- 
vey (White et al. Il997l ). but at significantly lower res- 
olution and signal-to-noise. However, the 1.4 GHz con- 
tinuum flux of 897 ± 148 /xjy beam~^ derived from the 
VLA -FIRS T data is significantly lower that that of Ivi- 
son (j2006[ ). To obtain a coherent dat aset, w e thus re- 
analyzed the A- array data by Ivison (|2006f ). We de- 
rive a continuum peak flux density of 1160 ± 33/xJy 
beam~^ (see Fig. 21 top left). This value is consistent 
with that found from the VLA-FIRST survey data, but 
significantly lower than that derived by Ivison (2006). 
Our result was independently confirmed by Ivison (2007, 
priv. comm.). We thus use our revised 1.4 GHz contin- 
uum flux in the following. 

3.5. Continuum Emission at 4-5 GHz 

We have detected continuum emission at 4.5 GHz to- 
wards APM 08279-^5255 using the VLA in B array. 
APM08279-f5255 remains unresolved in our observa- 
tions (see Fig. 21 bottom left). We derive a contin- 
uum peak flux density of 551 ± 40 /.t Jy beam~^. We 
also detect another continuum source about 3" south 
of APM 08279-f 5255, which is detected at a peak flux 
density of 180 ± 40/xJy beam~^ (this sourc e was first 
detected at 8.4 GHz by Ibata et al. Il99i^ see peak 
position derived below). At the given signal-to- noise, 
it is unclear whether this source is extended or not. 
APM 082794 -5255 has been detected at 4.5 GHz previ- 
ously (Ivison [2006| ). but at a factor of 4 worse linear res- 
olution (VLA C array, ~ 5.5"). Thus, APM 08279-1-5255 
and the source detected 3" south are blended at this 
resolution in Ivison's data, and the continuum flux of 
APM 08279-^5255 at 4.5 GHz is overestimated. The 
4.5 GHz continuum flux derived from Ivison's data (his 
Fig. 2, bottom panel) is in good agreement with the sum 
of the 4.5 GHz continuum fluxes of both sources in our 
analysis. We thus use our revised 4.5 GHz continuum 
flux of APM 08279-1-5255 in the following. 
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Fig. 4. — VLA maps of 1.4, 4.5, 8.4, and 23.5 GHz continuum emission towards APM 08279+5255 (1.4 GHz: robust 0, rest: natural 
weighting). The plus signs indicate the same positions as in Fig. [T] Top left: Image of the 1.4 GHz continuum emission at a resolution of 
1.42 "x 1.15" (as indicated in the bottom left corner). Contours are shown at (-3, 3, 5, 7, 9, 11, 15, 19, 23, 27, 31, 35)X(T (Ict = 33AtJy 
beam~^). This continuum map is integrated over a bandwidth of 43.75 MHz. Top right: Image of the 23.5 GHz continuum emission 
tapered to a resolution of 1.19"xl.l0". Contours are shown in steps of twice the rms noise, starting at ±3(t (Icr = 19 /ijy beam~^; 
19(T peak). This continuum map is integrated over a bandwidth of 100 MHz. Bottom left: Image of the 4.5 GHz continuum emission at 
a resolution of 1.55"xl.42". Contours are shown in steps of the rms noise, starting at ±2(t (Icr = 40 ^Jy beam"'^; 13(t peak). This 
continuum map is integrated over a bandwidth of 12.11 MHz. Bottom right: Image of the 8.4 GHz continuum emission at a resolution of 
0.73 " xO.42 ". Contours are shown at (-3, 3, 5, 7, 11, 15, 19, 23, 27, 31) Xc (Icr = 10 fi3y beam~^ ). This continuum map is integrated over a 
bandwidth of 100 MHz. In the 4.5 and 8.4 GHz maps, a radio source is detected about 3" south of APM 08279+5255 (at a = 08'^31™41=.576, 
5 = -|-52°45'14".66; 4cr and 7cr peak strength). This source remains undetected at 1.4 and 23.5 GHz (also undetected at 14.9 and 46.9 GHz). 



TABLE 3 
Spatially Resolved Components. 



Component 


Position 


5^ (CO) 


5^(8.4 GHz) 




(a;5) 




[MJy] 



APM08279(A) (08'^31™41=.708 ± 0=.001;+52°45'17".58 ± 0".02) 131 ± 16 205 ± 20 
(B) (08'^31™41^688 ± 0=.004;-|-52°45'17".34 ± 0".05) 87 + 16 204 + 20 
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Fig. 5. — Uniformly weighted 8.4 GHz continuum map. The 
beam size (0.63 "xO.22 ") is shown in the bottom left corner. Con- 
tours are shown in steps of the rms noise, starting at ±2(7 (Icr = 
20 /ijy beam~^). The source is clearly resolved into two peaks at 
this resolution (lOcr each), which coincide with the CO(J=l— »0) 
peaks (as indicated by the plus signs, see Fig.[T]l. 



3.6. Continuum Emission at 8.4 GHz 

Continuum emission from APM 08279+5255 was first 
detected at 8.4 GHz by Ibata et al. (jl99a) using the VLA 
in BnA array. We here combine this data with shorter 
B array spacings to improve the sensitivity of the data 
(see Fig. II bottom right). For APM 08279+5255, we de- 
rive a 8.4 GHz peak flux density of 446 ± 20 /iJy beam^^ 
(map convolved with a circular beam of 1.2" diameter 
for flux extraction), consistent with the value given by 
Ibata et al. (|1999t ) of 0.45 ± 0.03 mJy. We also detect 
the 'southern' source described in the previous subsec- 
tion. It is detected at higher signal-to-noise and res- 
olution in the unconvolved map than at 4.5 GHz. We 
thus use the 8.4 GHz data to derive the peak position 
of the source. The position of a=08^31™41^576±0^006, 
5=+52°45'14".66±0".03 is in good agreement with the 
peak of the 4.5 GHz continuum (source designation from 
now on: VLA J083141+524514). From the uncon- 
volved, naturally weighted emission line map, we derive 
a peak flux density of 83 ± 10/xJy beam~^. Note that 
VLA J083141+524514 is detected in none of the other 
bands, 3a upper limits are given in Tab.|4l 

We also imaged APM 08279+5255 with uniform 
weighting (Fig. [5]). The emission is clearly resolved into 
two peaks at this higher resolution, which coincide with 
the C0( J=l— s-0) peaks as indicated by the plus signs. In- 
terestingly, both peaks appear to have the same bright- 
ness at 8.4 GHz (205 ± 20/xJy beam~^ for the north- 
eastern peak, and 204 ± 20/zJy beam~^ for the south- 
western peak; corresponding to rest-frame peak bright- 
ness temperatures of 126±12 and 125+12 K). This is 
consistent with the naturally weighted image where the 
emission peaks centrally between the CO(J=1^0) peak 
positions (as both components are blended due to lower 
resolution). This different image brightness ratio relative 
to CO may be due to two non-thermal radio continuum 
components within the source with different spectral in- 
dices that are lensed differentially. VLA J083141+524514 
remains unresolved at this higher resolution. 
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Fig. 6. — GBT spectrum of C O(J=1^0) emission from 
APM 08279+5255 (Riechers et al.l2006D at a resolution of 5.86 MHz 
(75 kms~^, dashed line). The rms per channel is 65 ^tJy. The thin 
black line shows a Gaussian fit to the data. Overplotted is the 
peak flux densities of our new VLA CO( J=l— >0) data in the four 
120 kms~^ channels shown in Fig. |2] but tapered to the resolu- 
tion shown in Fig. \3\ The error bars indicate the rms per channel 
of 40 fi3y. This overlay shows that the GBT (single-dish) and the 
VLA (interferometer) detect the same amount of flux within the 
errors. 

3.7. Continuum Emission at 14-.9 GHz 

Continuum emission at 14.9 GHz was searched for to- 
wards APM 08279+5255. We tentatively detect emis- 
sion at the center of the CO(J=1^0) peak positions. 
This emission appears to be resolved and elongated along 
the same axis as the C0(J=1— >0) emission. The peak 
flux density measured at this position is 161 ± 44/iJy 
beam~^. However, since the detected flux density is <4ct, 
this measurement has to be treated with caution, and 
more sensitive data are required to confirm this result. 
Convolving this image with a circular beam of 1.2" di- 
ameter as above, we derive a flux density of 303+93 /iJy 
bcam^^, which we will treat as a tentative detection in 
the following. We summarize all our results in Tab. [2l 

4. ANALYSIS 

4.1. Spectral Energy Distribution 

Detailed continuum observations are now available for 
APM 08279+5255. Figure [5] shows the near-infrared-to- 
radio ('^1 /im-1 m, or 300 THz-300 MHz) spectral energy 
distribution (SED). The mid-IR-to-far-IR part of the 
spect rum has been model ed previously (Rowan-Robinson 
I2OOOI: Beelen et al. I2006D . indicating that a two compo- 
nent fit to the dust emission may be necessary to re- 
produce the data (a 'warm' AGN dust torus component 
plus a 'cold' starburst component : see a lso the multiple- 
component models by Blain et al. |2003[ ). Recently, Weifi 
et al. l|2007l ) have presented a model for the dust con- 
tinuum that also assumes two components, but which is 
additionally constrained by model parameters which re- 
sult from a study of multiple rotational transitions of CO 
and HON. 

However, no detailed analysis of the radio continuum 
properties of APM 08279+5255 has been conducted so 
far. We here present two modeling approaches. The first 
model aims at finding an overall good fit to the data 
while minimizing the number of free parameters. The 
second model aims at constraining which fraction of the 
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Fig. 7.— Redshifted CN(Af=l-+0) spectrum of 

APM 08279+5255 at a resolution of 3.5MHz (45 kms"!), 
obtained with the GBT. The velocity scale is relative to the 
redshifted CO(J=1->0) peak (2=3.911). The global rms is 65 /ijy. 
No emission is detected. The dashed lines indicate the range 
where the strongest fine structure components of the CN(7V=1— >0, 
^=0, J=3/2^1/2) multiplet, i.e. the F=3/2^1/2, F=5/2^3/2, 
F=l/2— >l/2, and _F=3/2^3/2 transitions, are expected to peak 
at 2=3.911. If the CN emission comes from a similar region than 
the CO emission, these components are expected to be blended. 

total radio continuum emission is due to star formation 
(i.e., not AGN-related"). 

4.1.1. Best Fit SED Model 

As our focus lies on fitting the radio data, we restrict 
our analysis to the mid-IR-to-radio part of the spectral 
data (60 /xm to 90 cm). For comparison only, we also 
show upscaled versions of quasar templates by Elvis et al. 
(fT994i) and Richards et al. (2006) on the optical to mid-IR 
part of the spectrum, which however are not particularly 
good fits to APM08279+5255's SED. 

Thermal dust emission is mainly caused by ionizing 
and non-ionizing ultra-violet (UV) radiation from mas- 
sive stars and (in quasars) possibly the central AGN. 
This radiation is processed by the dust and then re- 
radiated in the far-infrared wavelength regime. To min- 
imize the number of free parameters, we adopt the re- 
sults of Weifi et al. (2007^ to describe the dust contin- 
uum part of the SED. As discussed by WeiB et al., this 
imphes dust masses of Md.coid = 2.6 x 10^ ^^^M© and 

Ma warm 

— 7.0 X 10' /x^ M0 and dust temperatures of 
Td.coid = 75K and Td,warm = 220 K for the 'cold' and 
'warm' dust components, a dust disk with a magnified 
equivalent radius of vq = 1300 pc, relative area filling fac- 
tors of 78% and 22%, and a dust absorption coefficient 
exponent of /3=2 to describe its frequency dependence 
{fiL denotes the lensing magnification, the index d de- 
notes thermal dust emission). Using these parameters 
and equations 3, 4, 5, and 6 of WeiB et al. (2007) fully 
describes the model used here for the dust emission. The 
observed flux density of the dust emission will be denoted 
as S^^d in the following. To extend this model to the ra- 
dio continuum, we added a double power-law component, 
which we fit by the minimization technique. The full 
expression for the model flux density takes the form: 

Note that our models do not account for possible AGN vari- 
ability in this source (for which no evidence has been found to 
date). 



TABLE 4 

Continuum Fluxes of 
VLA J083141+524514. 



Frequency 




Configuration 








1.4 GHz 


<100 


A 


4.5 GHz 


180±40 


B 


8.4 GHz 


83±10 


BnA/B 


14.9 GHz 


<130 


B 


23.5 GHz 


<40 


BCD 


46.9 GHz 


<200 


CD 


Note. 


— Source 


is detected at 


position a 




Ml^576, S ^ 


-|-52°45'14" 


66. 





S.ii^) ^ S.^iiiy) + (1) 

2=fiat, steep 

where v indicates the frequency in GHz, c indicates 
a constant coefficient, a denotes a spectral index, and i 
indicates the index denoting the 'flat' and 'steep' power- 
law components. 

To solve equation (1), we attempted to fit the ob- 
served SED with three degrees of complexity. In the first 
step, we fixe d the dust term to the model parameters by 
Weifi et al. ()2007l see above), and fitted the 'flat' and 
'steep' components with a four-parameter double power 
law least-squares fit. The best fits to the 'flat' compo- 
nent result in asat values close to -0.1. However, the 
fitting uncertainty of this component is quite large due 
to the relatively small region where this component actu- 
ally dominates the emission. In the second step, we thus 
fixed aflat to -0.1. The remaining three free parameters 
are fairly well constrained. Our best-fitting model to the 
composite spectrum then gives a 'steep' power law index 
of astcop = —1.31 ±0.27, and 'steep' and 'fiat' coefficients 
of csteep = 1.09 ±0.13 and Cflat = 0.47 ±0.06. The 'steep' 
power law index suggests that APM 08279+5255 may be 
an ul tra-st eep spectrum source (asteep< —1-2; e.g., Bondi 
et al. l2007f ). which may indi cate t hat it evolves in a dense 
environment (Klamer et al. l20()6l ). In the third step, we 
left single or pairs of parameters of the Weifi et al. (200'$ 
dust model as free parameters to test the stability of the 
solution. This however did not further improve our best 
solution. Our final SED model is shown in the top panel 
of Fig. M 

4.1.2. Model SED of a Star-forming Galaxy 

Due to the fact that APM 08279+5255 is a compos- 
ite AGN-starburst system, it is likely that both com- 
ponents contribute to the continuum emission at radio 
wavelengths. Note that the 'warm' dust component in 
the Weifi etal. model of APM 08279+5255, which may 
be dominated by AGN heating, contributes ^90% to the 
far-infrared (FIR) luminosity. This finding implies that 
a clear separation of AGN and starburst contributions to 
the radio luminosity may be possible. 

We here investigate what fraction of the radio emis- 
sion can be produced by the starburst component in 
APM 08279+5255. Based on t he SE D of the Milky 
Way and first principles, Condon (1992) has developed a 
model for the radio-through-FIR emission from 'normal' 
galaxies, which is solely attributed to a stellar origin. 
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Fig. 8. — Near-infrared-to-radio spectral energy distribution (SED) of APM 08279+5255. The crosses indicate the measurements 
described in this paper, and the open square gives an estimate for t he 46.9 GHz cojitinuum flux. AU ot her signs indicate detections and 
upper limits from the literature (White et al. 119971 : Irwin et al. 119981 : Lewis e t al. 1199^ : Ibata et al. '1999'; Downcs ct al. 1999; Egami et al. 
[2000; Papadopoulos et al. 2001; Barvainis & Ivison 1200^ : Lewis et al. I2002al : Wagg etal. [2005, 2006; Boelon ct al. 2006; Ivison 2006; and 
Weifi et al. 2007). Top: The solid line represents a model fit to the thermal dust continuum, and a two-component power-law fit to the radio 
continuum emission. The components are indicated by the dashed lines. The parameters for the model for the dust part of the continuum 
are adopted from Weifi et al. C2007, ): the dotted lines indicate the 'cold' and 'warm' subcomponents of their model. The parameters of 
the 'fiat' and 'steep' spectrum components are determined by least-squares power law fitting. For comparison, upscaled versions of the 
'radio-quiet' quasar template (red) from Elvis et al. (1994) and the 'all' quasar template (blue) from Richards et al. (2006) are overplotted 
on the high-i' part of the SED. Bottom: The solid line represents a model of the thermal dust, free-free, and synchrotron emission SED 
components of a star-forming galaxy, as described by Condon (1992). The thick dashed lines indicate those components. Only the 'cold' 
dust subcomponent which is dominated by heating from star formation is considered for the model. The thin dashed line shows the model 
in the top panel for comparison (see text for more details). 
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This model contains three components: one for thermal 
dust emission, one for thermal (mostly free-free) radio 
emission, and one for non-thermal (mostly synchrotron) 
radio emission. The two radio components are described 
by power laws. This model also fits the starburst galaxy 
M82 very well. 

Free-free emission is either associated with the ionized 
gas in Hii regions surrounding hot young stars, or the en- 
vironment of the active galactic nucleus, and therefore an 
indicator for the total photoionization rate. If the ther- 
mal, optically thin free-free (or thermal bremsstrahlung) 
radiation emerges from a 10* K photoionized gas, a spec- 
tral index of ag = —0.1 (notation analogue to above) is 
expected, and thus assumed in the Condon model. The 
free-free power-law coefficient eg then parameterizes the 
emission measure, the source size and the electron tem- 
perature. However, note that free-free emission can only 
be described by a simple power law if the emission is 
optically thin. 

The (non-thermal) synchrotron emission is optically 
thin at typical radio wavelengths, and produced in su- 
pernova explosions and remnants (SNRs) through accel- 
eration of cosmic ray electrons, spiraling in a galaxy's 
magnetic field. Detection of strong synchrotron emis- 
sion from relativistic electrons thus indicates th e pres - 
ence of a large-scale magnetic field (Turner et al. Il998l ). 
In starburst galaxies, it may indicate a concentration of 
SNRs associated with localized and intense star forma- 
tion. Synchrotron emission usually shows a spectral in- 
dex of —1.2 < async < —0.4, with a typical value of 
-0.8, which is assumed in the Condon model. Here, the 
synchrotron power-law coefficient Cgync parameterizes the 
type II supernova rate (assuming that supernova rem- 
nants dominate th e non- thermal synchrotron emission). 

In the Condon l)1992[ ) model, all three components 
(dust heating, free-free, and synchrotron) are directly 
proportional to the star-formation rate (SFR) in a 
galaxy. The dust emissivity (3, the dust temperature Td, 
and the SFR are the only model parameters. We here 
adopt P=1.5 from Condon's model. Due to the fact that 
the 'warm' dust component in APM 08279-f 5255 may be 
dominantly heated by the AGN, we here only consider 
the 'cold' component, which likely traces the dominant 
fraction of the dust heated by star formation. We thus 
adopt T(3=75K. The only remaining free parameter thus 
is the SFR. Assuming a SFR of 4000 ^i^^ M© yr"! pro- 
vides a reasonable fit to the spectrum of the cold dust 
of APM 08279-1-5255 (see Fig. El bottom). However, this 
model onl y rec overs ~20% of the radio emission. Yun 
& CariUi pOO^ suggested some changes to the Condon 
model^'*, which would result in about a factor of 2 less 
radio fiux compared to the original Condon model. As a 
next step, we investigated if including the rd=220 K com- 
ponent (under the assumption that it is also heated by 
star formation) improves the fit. In this model, the peak 
of the dust spectrum limits the SFR from the 'warm' 
component to ~1500/i£^ Mq yr~^. Assuming such high 
values however requires to reduce the SFR of the 'cold' 
component to only few hundreds hJ^^Mq yr~^ in order to 
still reproduce the low-;/ end of the dust bump. This sig- 

^'^ They suggest a different initial mass function, an update to the 
Hcf normaUzation, and a small correction factor for the supernova 
rate. 



nificantly decreases the model-predicted radio flux, and 
thus gives an overall worse fit. The full radio flux could 
be recovered by including a very cold, ~10-15K dust 
component that is maintained by an extreme 15,000- 
20,000 /i^"'^ Mq yr~^ starburst. Such a component could 
exist in theory, but would require to be spread out over 
large scales in order to exhibit such low dust tempera- 
tures, and thus is basically ruled out by high spatial res- 
olution studies of this source. Although it is difficult to 
assess the validity of a simple Condon model for a warm, 
starbursting quasar hke APM 08279-1-5255, this result 
leads us to conclude that the radio continuum emission in 
this source may be dominanted by the AGN. This would 
be in agreement with our best-fit model, as it implies 
an ultra-steep radio spectrum component, which is usu- 
ally associated with a source where AGN-related emis- 
sion dominates the energy output at radio frequencies 
(e.g., high redshift radio galaxies and quasars; Athreya 
et al. Il997l : Pentericci etal. [2000). The ultra-steep spec- 
trum component in AGN-driven radio sources is often 
associated with extended radio emission. Interestingly, 
the best-fit model also requires a relatively bright 'flat'- 
spectrum component to originate from the radio AGN. 
The core of the archetypical radio galaxy Cygnus A has a 
flat spectrum at radio wavelengths, which dominates over 
the steep-spectrum emission from t he hotspots at higher 
frequencies (e.g., Eales et al. Il989h . A similar 'compos- 
ite' effect (extended steep-spectrum component, flat core 
component) may cause the apparent flattening of the ra- 
dio spectrum of APM 08279+5255. 

4.2. Comparison to Previous CO(l-O) Imaging 

The CO(J=1^0) transition line in APM 08279+5255 
was targeted twice before with the VLA, first at Car- 
ray resolution (1.5"xl.4", Papadopoulos et al. l2001h . 
and then followed up at higher B array resolution 
(0.39"x0.28", Lewis et al. 12002a'). Both observations 
were setup in quasi-continuum mode, i.e. lacking any 
spectral information of the C0(J=1— >0) line. In both 
observations, C0(J=1— >0) line and 23.3649 GHz contin- 
uum emission were detected towards the source. Pa- 
padopoulos et al. (2001) detect unresolved continuum 
emissi on at a flux density level of 0.30 ± 0.05 mJy. Lewis 
et al. (|2002af ) detect continuum emission at a fiux density 
of 0.41 ± 0.07 mJy, and the resolved structure coincides 
with our continuum map in structure and position. Both 
flux densities are in agreement with our result within the 
errors, although some intrinsic variability cannot be ex- 
cluded. 

Papadopoulos et al. (|2001l ) report the detection of re- 
solved C0( J=1^0) emission towards APM 08279+5255, 
which extends over scales of >7"x2.25" (their 3cr con- 
tour) at a resolution convolved to 2.25 "x 2.25" (7" cor- 
respond to ~30kpc at the target redshift of 3.911). They 
quote an integrated fiux for the 'nuclear CO(J=1-^0) 
emission' (the inner ~1") of 0.150 ± 0.045 Jy kms"^ 

From our new VLA observations, we derive an inte- 
grated CO(J=1-^0) flux of 0.168 ± 0.015 Jy kms'^ as- 
suming a line FWHM of 556 ± 55 kms^^ (Riechers et 
al. [20061: width of our bandpass: 558 kms ^). This cor- 
responds to a CO(J=1^0) line luminosity of ^^0(1-0) 
= (10.6±0.9)xl0iVZ^K kms-ipc^ (see Tab.©. From 
the CO(J=1-^0) spectrum obtained with the GBT, 
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Riechers et al. ()2006D find 0.152 ± 0.020 Jy knis"^ The 
3cr contour of the extended C0( J=l-^0) reservoir of Pa- 
padopoulos et al. (|2001l . their Fig. 1) corresponds to a 
flux density of 120/xJy beam~^. In our CO(J=1^0) 
map which was tapered to a linear resolution of ^1.2" 
(Fig. [HI similar to the resolution of their observations), 
this would correspond to 6.3a. Even at less than half this 
flux level (3CT=57/iJy beam~^), we find no evidence for 
any extended flux. We note that our observations even 
incl ude sh orter (D array) baselines than Papadopoulos et 
al. (IMl and are therefore more sensitive to extended 
structure. We thus exclude the possibility that our obser- 
vations outresolve the extended emission, and conclude 
that there is no bright, extended C0(J=1— >0) reservoir 
in APM 08279+5255. 

Lewis et al. (,2002a, ) find resolved CO(J=1^0) emis- 
sion on sub-arcsec scale at moderate signal-to-noise. The 
peak of the emission reported in their paper roughly co- 
incides with our north-eastern C0(J=1^Q) peak. How- 
ever, our new, more sensitive (by a factor of 2 in terms 
of rms noise) high-resolution maps do not show the ex- 
tended structure that is present in their C0(J=1— >0) 
map, which we thu s conclude to be a noise artifact. 
Lewis et al. (jlOOl^) derive an integrated CO(J=1^0) 
flux of 0.22 ± 0.05 Jy kms"-'^ from their observations, 
in agreement with the three independent results given 
above within the errors (adopting only the 'nuclear 
C0(J=1— *0) emission' from the report of Papadopou- 
los et al.). 

4.3. On the 4^6.9 GHz Continuum Emission 

As mentioned previously, the continuum emission at 
the line frequency (46.9 GHz) was not monitored in paral- 
lel during the (ar chival ) C0(J=2— >1) observations. Pa- 
padopoulos et al. (2 0011 ) report a limit of 3<t = 0.9 mJy on 
the continuum emission at 43.3 GHz. This corresponds 
to ~55% of the C0( J=2^1) peak flux. Therefore, the 
'real' C0(J=2^1) flux towards APM 08279+5255 may 
be lower by up to a factor of 2. Based on our SED 
modeling, we now attempt to set additional limits to 
the continuum contribution to the emission detected at 
the C0(J=2^1) frequency to attain some tighter con- 
straints. 

A first-order estimate for the 46.9 GHz continuum flux 
is provided by examining the overall spectral energy dis- 
tribution of APM 08279+5255. The SED model pre- 
sented in Fig.[8]predicts a continuum flux at 46.9 GHz of 
405 /iJy. This estimate is plotted as the square symbol 
in Fig. [S] A lower constraint is provided by the fit to 
the thermal part of the SED (75/iJy), and an upper con- 
straint is given by linear interpolation of the 23 GHz and 
94 GHz continuum flux measurements (715 /iJy). For il- 
lustration, these constraints are plotted as error bars to 
this model-predicted data point. We however note that 
the 23 GHz and 94 GHz measurements were obtained at 
different epochs, therefore this estimate does not account 
for possible variability of the millimeter continuum. 

As another estimate, we assume that both the emission 
from CO(J==1-^0) and C0(J=2^1) are fully thermal- 
ized and optically thick. The C0(J=2— >1) peak flux 
should then be four times higher than the C0(J=1— »0) 
peak flux. For the C0(J=1— s-0) peak flux, an aver- 
age over the whole 558 knis^^ bandpass (210+19 /iJy) 
is assumed to match the velocity resolution of the 
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Fig. 9. — Model spectrum of the CN(A''=1^0) emission in 
APM 08279+5255. The horizontal axes indicate the rest and ob- 
served frequencies. The vertical axis indicates the predicted inten- 
sity, normalized to the peak of the emission. The relative intensi- 
ties of the fs and hfs components are computed for the LTE case. 
The crosses indicate the hfs components. The dashed lines are 
Gaussian fits to the summed intensities in the two main frequency 
bins, assum ing the average width of the higher J CO lines (WeiB 
et al. 120071 ) for the subcomponents (labeled '1' and '2'). The solid 
line is a sum of all components, and indicates the model-predicted 
CN(Af=1^0) Une shape. 

C0( J=2^1) observations as closely as possible. In com- 
parison to the measured flux at the C0(J=2^1) fre- 
quency (1637+65 /iJy), this predicts a contribution of 
the continuum flux of ^800 /iJy, which is twice as high 
as the SED-based estimate. This discrepancy may be ex- 
plained by different optical depths of the emission from 
different rotational CO transitions. If the opacity of the 
C0(J=1— >0) emission is much lower than that of the 
higher J transitions, it may actually have a lower lu- 
minosity than the emission from these transitions [con- 
trary to the effects of (sub-)thermal excitation]. This 
would allow the C0(J=2— >1) line peak flux to be more 
than four times higher than the CO(J=1^0) peak flux. 
Note that such an effect may actually be observed to- 
wards APM 08279+5255 in the J>4 transitions, as the 
CO( J=4-> 3)/CO(J=1^0) luminosity ratio is >1 (Weifi 
et al. I2007t see also Tab. O. 

We conclude that APM 08279+5255 likely exhibits a 
46.9 GHz continuum flux of at least 0.4 mJy, i.e. 25% 
of the measured C0(J=2^1) flux. We thus estimate 
the 'real' C0(J=2— >1) line peak flux to be ~1.2mJy. 
From this, we derive a C0(J=2^1) line luminosity of 

^co(2-i) = (12-8±2.9)xl0iVZ^ K kms'^ pc^. 

4.4. A Limit for the CN Luminosity 

Due to its fine structure (fs) and hyperfine structure 
(hfs) splitting, the CN(iV=1^0) emission line is dis- 
tributed over 9 lines^^ in its v=Q vibrational state. These 
lines are distributed in 2 main components separated by 
about 350 MHz in the rest frame. While the hfs transi- 
tions within these main components are blended if their 

We assume that CN is in its ground electronic state (^S), and 
that the spins couple according to Hund's case (b) coupling scheme: 
N + S = J, {is coupling) and J + Iq + = F (hfs coupling). Here, 
N is the rotational angular momentum vector, S is the electronic 
spin, and / is a nuclear spin. 
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TABLE 5 

Line Intensities, Luminosities, and Ratios. 



Line 




/ 


L' 


-'^'/-^'co(l-O) 


data reference 






[Jy kms"'^] 


[10^° K kms^i pc2] 






CO{J=l- 


>0) 


0.168 ± 0.015 


10.6 ± 0.9 


1 


1,2 


CO(J=2- 


*1) 


0.81 ± 0.18 


12.8 ± 2.9 


1.21 ± 0.30 


1,3 


CO(J=4- 


»3) 


3.7 ± 0.2 


14.7 ± 0.9 


1.39 ± 0.16 


4 


CO(J=6- 


»5) 


6.7 ± 1.2 


11.8 ± 2.1 


1.12 ± 0.23 


4 


CO{J=9- 


>8) 


11.8 ± 0.6 


9.2 ± 0.4 


0.89 ± 0.10 


4 


CO(J=10- 


-9) 


11.9 ± 2.0 


7.5 ± 1.2 


0.71 ± 0.14 


4 


CO(J=ll- 


--10) 


11.3 ± 1.9 


5.9 ± 1.0 


0.56 ± 0.11 


4 


CN(Af=l- 


-0) 


<0.06 


<3.9 


<0.37 


1 



Refer ences. — [1] this work; [2] Rieciicrs ct al. l2006t [3] Papadopoulos ct al. 120011 : [4] WciB ct 
al. [20071 . 

Note. — Luminosities are apparent values uncorrected for gravitational magnification. They are 
derived as described by Solomon et al. I I1992D : i' [K km s~^pc^] = 3.25 X lo'' Iv~^^DI{1 + z)~^ fj.^^ , 
where I is the velocity-integrated line fiux in Jy kms~^, Dl is the luminosity distance in Mpc, 
fohs is the observed frequency in GHz, and /^i is the lensing magnification factor. 



kinematical broadening is similar to that of the CO /HCN 
lines in APM 08279-1-5255, the components themselves 
are still clearly separated. To obtain a proper limit on 
the integrated CN(A^=1^0) line flux (and thus on the 
total CN luminosity) , we calculated a synthetic line pro- 
file, assuming optically thin emission in the Local Ther- 
modynamic Equilibrium (LTE) case to derive the relative 
intensities of the hfs components. The relative intensi- 
ties of the components were c ompu ted using the labora- 
tory data from Skatrud et al. (|1983( ). and ap proximation 
(4) of equation (1) of Pickett et al. ()1998[ ). Assuming 
that all CN(A^=1^0) hfs components are kinematically 
broadened in the same way as the high er J CO lines 
(-500 kms-i FWHM; see Weifi et al. 1200% . we ob- 
tain the synthetic line profile displayed in Fig. [S] (solid 
line). The two main components (dashed lines) have a 
peak strength ratio of 2:1. The brightest main compo- 
nent is fully covered by our bandpass, the second com- 
ponent, however, was close to the edge, and thus only 
covered in part by the part of the spectrum that was 
properly calibratable. For the derivation of the limit on 
the integrated CN(iV=l— >0) line flux, we thus consider 
the peak flux limit for the brightest component to derive 
the relative intensity limits for all subcomponents in the 
LTE case. Under these assumptions, we derive a limit 
of /cN<0.06 Jy kms~^ for the integrated line flux, and 
thus a limit of Lcj,j<3.9x 10^° K kms~^pc^ for the 
line luminosity. This corresponds to a luminosity ratio 
limit of L'(.^/L'^Q<0.37. Note that Guelin et al. (21)0?) 
recently observed the CN(A^=4— >3) transition towards 
APM 08279-1-5255; however, due to strong blending with 
the HNC(J=5^4) line, it remains unclear whether or 
not the line is detected. 



corrected), about 10 times the radio luminosity of M87 
(-^i.4 GHz = (1.76±0.07) X 10^^ W Hz-^ Laing & Peacock 
Il980f ). In an analogous manner, we derive a 5 GHz radio 
luminosity of i5GHz = (4.75±0.34)xl025/iZ^ WHz-i. 
Here, the model prediction would correspond to a 
spectral index of only a5=-0.67. 

Studies of the radio luminosity function of quasars 
and its evolution have identified a bimodality in the 
source distribution, separating them into radio-quiet and 
radio-loud quasars. However, two generally different def- 
initions of radio loudness have been put forward, one 
based on the optical-to-radio flux ratio (e.g., Schmidt 
119701 : Kellerman et al. ,1989 ) , and one only based on 
the m onochromatic r adio l uminosity (e.g.. Peacoc k et al . 
119861: Schneider et al. ll992|) . However, Ivezic et al. (|2002D 
argue that for optically selected quasars, both defini- 
tions are similar within the flux-limited samples that 
were examined. We thus restrict our analysis to the 
radio luminosity-based estimate here. Due to different 
assumptions and definitions, a range of threshold val- 
ues for the radio loudness definition is found in the lit- 
erature. For the radio-based definition, we here adopt 
the range of L5GHz^1-0-5.6)x10^'^ WHz~^ as given by 
Schneider etal. (Il992f ). Changing the cutoff within this 



range does not change the relations they find dramati- 
cally. Even without correcting for gravitational lensing, 
APM 08279-f 5255 clearly falls within the transition re- 
gion between radio-quiet and radio- loud sources. As this 
source is known to be substantially gravitationally mag- 
nified (see also discussion below), we conclude that it is 
a radio-quiet quasar by this definition. 



4.5. Radio Luminosity 
In general, the monochromatic rest-frame 1.4 GHz ra- 



dio luminosity is defined as: 



4 GHz 



(1 



U+ai.4 



^1. 



4 GHz- 



(2) 



However, we here calculate Li.4GHz based on the 
(rest- frame) 1.4 GHz flux predicted by our best-fit SED 
model. In the above equation, this would correspond to 
using a spectral index of ai.4=-1.05. We thus obtain 
-/-i.4GHz = (l-90±0.05)xl02VZ^ WHz-i (not lensing 



4.6. Morphology at X-Ray to Radio Wavelengths 

Optical imaging of APM 08279-^5255 with the Hubble 
Space Telescope (HST) Near Infrared Camera and Multi- 
Object Spectrometer (NICMOS) has revealed three 
pointlike, almos t coUi near images with very similar col- 
ors (Ibata et al. Il999f ). A spectroscopic follow-up study 
with the Space Telescope Imaging Spectrograph (STIS) 
spectrograph on board HST has shown that these three 
images have essentially identical spectral shapes, and 
thus are, indeed, all imag es caused by gravitational lens- 
ing (Lewis et al. l2002bf ). These three optical images 
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Fig. 10. — Morphology of APM 08279+5255 at different wavelengths. Left: Contours of 2.2 /im continuum emission (Egami et al. I2OO0I) 
overlaid on the CO(J=l— >0) emission. Middle: Contours of 814 nm continuum emission (HST archive) overlaid on the CO(J=l— >0) 
emission. The inset shows the optical image only. Right: Contours of hard X-ray emission (Chartas et al. 12003) overlaid on the CO( J=l— >0) 
emission. Contours in all panels are 20%, 30%, 90% of the peak flux of the emission, but with 23.3% and 26.7% added in the middle 
panel to emphasize the third image. The CO greyscale scales between 2a and the peak flux of the emission (see Fig. [I}. The crosses 
indicate the CO peak positions. The size of the synthesized beam of the CO observations is indicated in the bottom left corner. 



A, B, and C^^ are separated by d(A,B)=0.377"±0.002" 
and d(A,C)=0.150"±0.006", and have brightness ratios 
of &(A,B)=0.773 ± 0.007 and 6(A,C)=0.175 ± 0.008, 
i.e., &(A-t-C,B)~0.66. The north-eastern CO peak cor- 
responds to a blend of optical images A-l-C, and the 
south-western CO peak to optical image B. The CO im- 
age peak brightness ratio of 6(SW,NE)=0.66 ± 0.15 thus 
agrees very well with the optical observations. The sep- 
aration of the CO images is by ~20% lower than that of 
optical images A and B, which agrees with the assump- 
tion that the north-eastern CO component is actually a 
blend of optical images A-l-C, which we will assume in 
the following. 

In their analysis, Ibata et al. ( |1999| ) find an offset of 
0.6" between the positions of HST NICMOS image A 
and the 8.4 GHz continuum source, which is not thought 
to be real. From our analysis, we find that the obser- 
vations in six different radio frequency bands are con- 
sistent in position among each other to high precision 
(and with the radio positions from previous studies). We 
also find that the offset between the radio position and 
a more recent image taken with HST Wide-Field Plane- 
tary Camera (WFPC2; HST archive, see also Fig.[Tni) is 
<0.2". This imphes a large offset between the WFPC2 
and NICMOS images, which are close in wavelength and 
likely dominated by compact, parsec-scale AGN emis- 
sion which is expected to be spatially coincident. In the 
following, we thus assume that these offsets are due to 
HST astrometric calibration errors, and shift all images 
at other wavelengths to the radio coordinates. 

The three images of APM 08279-^5255 were also re- 
solved in the (observed frame) near-infrared (NIR) con- 
tinuum at 2.2 /im in ground-based adaptive opt ics ob ser- 
vations with the Keck telescopes (Egami et al. |2000[) at 
a resolution where images A and C are still significantly 
blended, similar to our C0(J=1— >0) observations. The 
left panel of Fig. [TUl shows an overlay of the 2.2 /im con- 
tinuum emission and the C0(J=1— >0) emission. The 
2.2 /xm has b een d e-rotated (see caption of Fig. 1 of 
Egami et al. |2000| for the rotation angle) and aligned 

By convention, A is the brightest image, and C is the faintest 



to the radio coordinates. The resolution of the NIR ob- 
servations is by about a factor of 2 higher than that of 
the CO observations, and thus peaks close to the po- 
sition of image A rather than between A and C. The 
image separation of the two peaks thus is d(A,B)=0.38", 
in agreement with the optical images (Fig. [TOl middle 
panel). Taking this into account, it is striking how sim- 
ilar the morphologies of the NIR and CO emission are, 
even though they are expected to emerge from largely 
different scales (hot dust around the central AGN vs. 
cold molecular gas in an extended star- forming ring), and 
thus to be differentially lensed. This similarity of mor- 
phologies even holds for hard X-ray emission from the 
central engine (Chartas et al. l2002l see Fig. [lOl right), 
which is expected to be by many orders of magnitude 
more compact than the CO emission. Even the image 
separation (i(A,B)=0.38"±0.01" of the two main peaks 
(image A and C are again significantly blended) is the 
same as in the optical/NIR, and thus in agreement with 
the CO observations under the above assumptions. 

In addition, the rest-frame 2.6 mm continuum under- 
lying the CO(J=1^0) line shows a different, more ex- 
tended structure than the CO emission^^ (see Fig. [T]). 
This may be due to a foreground source that emits in 
the continuum, but does not contribute to the line emis- 
sion due to its different redshift. Another, probably more 
likely explanation is that the different structure is due to 
differential lensing. If the 2.6 mm continuum emission 
were to be due to star formation, it would be possible 
that it is dominated by a free-free component that is not 
co-spatial with the CO emission, but extended out to kpc 
scales. On the other hand, the SED of APM 08279-f 5255 
indicates that a major fraction of emission at this wave- 
length may also be AGN-related. If correct, this would 
suggest that the differentially lensed emission is due to 
100 pc to kpc-scale outflow from the nucleus. Higher res- 
olution observations using a very long baseline radio in- 
terferometer are required to further investigate the na- 
ture of this phenomenon. Also, imaging the rest-frame 
FIR continuum at comparable resolution may shed more 

^"^ A possible caveat is that the bandwidth of our observations 
is too narrow to image the CO{J=l— >0) linewings. It however is 
highly unlikely that this significantly alters this conclusion. 
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light on this issue. In the following, we will refer to this 
emission as 'the extended component'. This differentially 
lensed, extended component may also contribute signif- 
icantly to the continuum emission at 7.3 mm (8.4 GHz 
observed frame), which also shows a different structure 
compared to other wavelengths. 

The morphological similarity of the strongly lensed 
quasar APM 08279-1-5255 on largely different scales 
(<0.1pc to few 100 pc) as well as the different struc- 
ture of the extended continuum component (probably 
few 100 pc to kpc) severely constrains the lens configu- 
ration, and any valid lens model will have to be able to 
reproduce this scaling effect. 

5. GRAVITATIONAL LENSING 

5.1. Direct Constraints on the Lensing Properties from 
the CO Observations 

The brightness temperature Tb of a lensed source is 
a conserved quantity under gravitational lensing, as the 
latter is a purely geometrical effect. Due to the fact that 
APM 08279-^5255 is resolved in our CO(J=1^0) obser- 
vations, we can derive the observed brightness tempera- 
tures T^^^^'^ at redshift z for images A-l-C and B directly 
from the observed peak fluxes using 

rpOhs,Z C , . 

where ilbcam is the beam solid angle. For Gaussian 
beam and source shapes, this equation can be simplified 
to 

T°^'^'- ~ 1.359 X 10' % S,, (4) 

where Tb is given in K, Aobs is given in mm, the beam 
major and minor axes ^maj and ^min are given in arcsec, 
and Siy is given in /iJy. These brightness temperatures 
can be converted to rest-frame observed brightness tem- 
peratures T^^^ via 

T^^' = {1 + z)T°^'''\ (5) 

While CO(J=1^0) images A-l-C and B are individu- 
ally detected in our observations, their true size is likely 
smaller than our observing beam. The observed bright- 
ness temperatures T^^^ thus are diluted by the beam 
(radius: rbeam) and smaller than the true Tb. The size 
of the lensed images can be expressed by an equivalent 
radius tq, which assumes that each image is a filled, cir- 
cular lensed disk. Energy conservation then gives: 

T^'^yn = (ro/rbeam)'. (6) 

This relation is implied in the plot shown in Fig. Illh . 
The observed, beam-diluted (Rayleigh- Jeans) brightness 
temperatures of images A-l-C — 23.2 K) and 

B [T°''^'" = 12.6 K) are indicated by vertical lines. The 
three curves indicate the beam-corrected brightness tem- 
peratures for different equivalent radii (which, by def- 
inition, cross with the observed rest-frame brightness 
temperatures at rp /rbeam j as indicated by the long hor- 
izontal line) for images A-fC, B, and all images to- 
gether. From Large Velocity Gradient (LVG) model- 
ing of the CO(J=1->0) to Cb( J^ll^lO) SLED, Weifi 
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Fig. 11. — Observational constraints on the intrinsic brightness 
temperature Tj, and the lensing properties of APM 08279+5255. 
a: Relation between Tj, and the lensed equivalent disk radius rg. 
The vertical lines indicate the observed CO(J=1^0) brightness 
temperatures for the blended images A and C and image B, and 
the LVG predicted brightness temperatures for the 'cold, dense' 
(50 K) and 'warm' (60 K) CO(J=l— ►0) components from the model 
by Weifi et al. (2007). The horizontal lines indicate the lensed 
equivalent disk radii for the different images for a true T^, of 50 K, 
and the total ro for all images. The fourth horizontal line indicates 
the VLA beam radius, and also appears in panels b and c. b: 
Relation between the true, unlensed equivalent disk radius r^^^'^ 
and the effective lensing magnification factor fi'j^ for different T^. 
For Tb=50K, different fi'^ are highlighted, c: Relation between 
the true, unlensed radius rtruc of the CO ring and ^"j^ for the 
T]-|=50K solution assuming an inner ring radius r^^ of 50 pc for 
different inclinations i and area filling factors /a. For comparison, 
a solution without an inner boundary is shown (see text for details). 
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et al. (|2007f) have derived that the molecular gas in 
APM 08279+5255 can be modeled with two gas com- 
ponents, a 'cold, dense' component with a H2 density 
of n(H2)=10^ cm~^ and a kinetic gas temperature of 
2kin=65K which contributes about 70% to the total 
C0( J=l-^0) luminosity, and a 'warm' component with a 
H2 density of n(H2)=10'* cm~'^ and a kinetic gas temper- 
ature of rkin=220K which contributes about 30% to the 
total CO(J=1^0) luminosity. For the 'cold, dense' com- 
ponent, their model predicts a true CO(J=1^0) bright- 
ness temperature of T™''^=50K, and for the 'warm' com- 
ponent, it predicts r™™=60K. Both Tb are indicated 
in Fig. Illb by vertical lines. Assuming rb=50K, as pre- 
dicted by the LVG models for the dominant gas compo- 
nent, images A+C and B fill equivalent disks with radii 
of r^'''^=717pc and Tq =530pc. This corresponds to a 
total lensed equivalent disk with a radiu s of rn=892 pc, 
consistent with the results of WeiB et al. ( 2007) within the 
uncertainties. This also implies beam area filling factors 
of 46% and 25% for images A-l-C and B. These are sig- 
nificant fractions of the beam size, and may indicate that 
increasing the resolution by only a factor of 2-3 would 
be sufficient to detect substructure within the individ- 
ual images. Assuming a 30% contribution of a warmer 
component with Tb=60K does not significantly alter any 
of these predictions. Even assuming much higher T^*-* 
than typically observed in the central regions of nearby 
warm, ultra- luminous infrared galaxies (ULIRGs, e.g., 
Downes & Solomon 1998) of a few hundred Kelvins does 
not change the predicted vq by more than about a factor 
of 2. Also, the directly observed (beam-diluted) bright- 
ness temperatures limit the true Tb to be at least ~30 K 
considering that the contribution of image C to A-l-C is 
likely minor. 

The effective lensing magnification is defined as 
the ratio between the observed apparent luminosity and 
the true luminosity. It thus is a direct measure for the 
ratio of the observed equivalent disk size tq and the true 
equivalent disk size r^'^^" via 

A^f = iro/r'jn'- (7) 

This relation is displayed in Fig. [Tib for a range of Tb. 
The Tb=50K model discussed above is indicated by the 
solid line and the boxes (highlighting selected magnifi- 
cation factors). Using only the given constraints with- 
out modeling the system in more detail, a large range of 
magnification factors would be in agreement with the CO 
data. However, the impact of changing the intrinsic Tb 
is relatively minor, as all curves occupy a narrow range 
in this plot. 

The case of a completely ffiled, circular CO disk as dis- 
cussed so far is a very particular solution, even appreciat- 
ing the fact that the CO is likely situated in a relatively 
compact, rotating circumnuclear ring. In a more realis- 
tic approach, the CO disk has an inner boundary rin at 
a certain distance from the central hot nucleus, due to 
the simple fact that molecular gas cannot survive at dis- 
tances where the ionizing radiation field is strong enough 
to dissociate the diatomic CO (and H2) molecules. It 
thus is a ring rather than a disk. Also, it is likely that we 
do not see this ring face-on, but at a certain inclination i 
towards the line of sight (LOS). Also, the molecular ma- 
terial may be distributed in clouds and/or clumps rather 



than smoothly, and thus may have an area filling factor 
/a lo wer th an 1. This can be expressed as (see also Weifi 
et al. [2007l) : 



Note that i—0° corresponds to face-on. Also, we use 
a thin disk approximation, neglecting the geometrical 
thickness of the CO ring. This simplification leads to an 
overprediction of rtiuo for high inclinations. This approx- 
imation is justified, as the high observed optical bright- 
ness of the source in combination with the detection of 
large amounts of dust in its circumnuclear torus suggests 
that we have a relatively unobscured view on the AGN, 
and thus that we likely see the galaxy at a relatively low 
inclination (i<35'^). 

In Fig. Illb . a range of inclinations is shown for the 
Tb=50K model in the nf^-rn^c plane. An inner ring 
boundary of rin=50pc is assumed unless stated other- 
wise. For comparison, a solution without an inner ring 
boundary is shown for i=0° . The change of the predicted 
rtruc is only minor, in particular for the solutions with 
low magnifications. In general, the figure shows that the 
range of solutions for high inclinations is significantly 
different from that for low inclinations. However, in the 
preferred range of 0° < i < 35°, the solutions are fairly 
similar. While rin and i can be constrained relatively 
well by existing observations, this is not at all true for 
/a. As an example, a solution with i=75° and /a=0.5 
is shown. The difference in area filling factor has a large 
impact on the predicted solutions. As another exam- 
ple, solutions for i—0°, /a=0.5 and i=60°, /a=1.0 are 
indistinguishable. The structure of the molecular ISM 
thus is an unknown in this discussion that has signif- 
icant impact. There are two more interesting results 
to note: First, the resolution of our observations and 
A+C/B brightness ratio alone set a meaningful lower 
limit on the magnification of models with /a cos i signif- 
icantly lower than 1. Second, the models with /acosi 
close to 1 predict very compact CO rings for high mag- 
nifications (rti.uc=80-125pc). 

Although the results presented in this subsection sig- 
nificantly constrain the allowed parameter space for po- 
tential lensing models, they do not allow to strongly con- 
strain the actual magnification. In the following, we thus 
discuss existing lensing models in more detail, show their 
limitations, and suggest a new model which overcomes 
part of the problems of previous models. 

5.2. Lens Modeling: Previous Models 

In most classical lensing models, the surface density of 
the lensing galaxy is expressed by an ellipsoidal with a 
core of finite radius. However, high-resolution observa- 
tions of the central regions of galaxies indicate that the 
luminosity profiles of some of the most massive galaxies 
appear to have central cu sps ra ther than cores with fi- 
nite radii (e.g., Faber et al. 119911) . Such distributions are 
reminiscient of the cusp power-law density profiles with 
a break radius used for dark matter halos in cosmological 
simulations (e.g., Navarro, Frenk & White 1997, 'NFW'), 
but, in these observations, are seen for the stellar com- 
ponent of the galaxies. The (baryonic) density profile of 
the lensing galaxy of APM 08279-1-5255 has been mod- 
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eled with both core (Ib ata at al. ll999t Eg ami et a l. l2000f) 
and cusp (Munoz et al. 120011 : Lewis et al. [20023) config- 
urations previously. 

5.2.1. Cored Singular Isothermal Ellipsoid Models 

The model ol Egami et al. ()2000[ ) is a modified singu- 
lar isothermal ellipsoid (SIE) model with a finite core. 
Such models have six parameters: the Einstein radius, 
the core radius, the ellipticity ol the lensing potential 
and its position angle, and the (2-dimensional) center 
position ol the source. For a three image source like 
APM 08279+5255, observations olFer six relevant con- 
straints for the model: two relative image brightnesses, 
and two relative (2-dimensional) image positions. In ad- 
dition, the fact that the three images are almost coUinear 
requires the ellipticity to be close to zero, and thus pre- 
dicts an almost circular lensing potential. The geom- 
etry of lensing systems requires that any non-singular 
mass distrib ution produces an odd number of lens im- 
ages (Burke [l981|) . However, if the core is very small, 
one of the images is de-magnified, and thus unlikely to 
be observable (note that an even number of images is ob- 
served toward most lens systems) . The fact that the third 
image of APM 08279+5255 is not strongly de-magnified 
thus also restricts the core to be non-singular, and thus 
to have a ra ther large radius within this model. Egami 
et al. (I200I find a core radius of 0.2", corresponding 
to 1.2-1.7kpc at 0.5<z<3.5, the likely redshift range of 
the lensing galaxy, which appears extreme even for the 
most massive elliptical galaxies. In particular, quasar 
image C is located only 0.03" (180-250 pc) from the cen- 
ter of the lensing potential. As the authors already state 
themselves, one thus might expect differential reddening 
of image C in this configuration, which is not observed. 
Howev er, the most problematic point of the Egami et al. 
(|2000f ) model is that it predicts a rapid change of mor- 
phology with source size (their Fig. 7). For the X-ray 
to near-infrared observations, where most of the emis- 
sion likely comes from the central <0.1-lpc region, their 
model correctly predicts a three image configuration. 
This even holds true for a source size of 20 pc, where the 
images are elongated, but still clearly seperated, and thus 
would be picked up as three pointlike images by obser- 
vations at a resolution of, e.g., 0.3". However, for source 
sizes of 50 pc and larger, the images start to form arcs 
and rings, and ultimately at 220 pc, a filled sphere. How- 
ever, in our C0(J=1— >0) observations, we probe scales 
of 100-500 pc (see above), depending on the lensing mag- 
nification, but still see the 3-image structure as observed 
in the X-ray-to-NIR wavelength regime. Moreover, the 
extended components seen in radio continuum emission, 
which fikely probes few 100 pc to kpc scales, gets lensed 
into an arc-like structure rather than a filled sphere (see 
Fig. [Tl right). While this model predicts high effective 
lensing magnifications of ^f^'^lOO (and thus a plausi- 
ble explanation for the extreme observed properties of 
APM 08279+5255), it is ruled out by the structure de- 
tected in our CO and radio continuum maps. 

5.2.2. Cusp Models 

As mentioned above, there is mounting evidence that 
galaxy cores are mostly very small, and that their mass 
distribution possesses a significant cusp (Rusin & Ma 



I2OOII: Winn et al. I2004D . Munoz et al. (|2001h thus sug- 
ges ted a model complementary to the one by Egami et 
al. (|2000|) . assuming a cusp configuration for the lens in 
APM 08279+5255. Instead of a core radius, such mod- 
els are parameterized by a cusp power-law index and a 
break radius, i.e., 7 parameters in total, and thus are 
underconstrained by a three-image lens without further 
assumptions. For APM 08279+5255, cusps with a small 
break radius are favoured, and lead to high magnification 
factors of /i^*^'--^100, similar to the Egami et al. (2000) 
model. However, the source can only be modeled with 
relatively shallow cusps with a small power-law index 
(0.2-0.4 rather than 1-2 as observed toward most galax- 
ies). Such shallow cusps are, indeed, relatively similar 
to cores with finite radii (corresponding to a power-law 
index of 0). The authors conclude that the lens may be 
a spiral galaxy, since some spiral galaxies appear to have 
finite central densities rather than steep cusps. 

Kee ton & Kochanek ()1998l ) and Bartelmann & Loeb 
|1998!) have shown that highly fiattened, highly inclined 
potentials, as found, e.g., for disks of edge-on spiral 
galaxies, can form a 'naked cusp' caustic. In such a con- 
figuration, the inner diamond caustic extends outside the 
elliptical caustic, producing three roughly coUinear im- 
ages with similar brightness. The central third image is 
not strongly de-magnified, which is different from typical 
core and cusp models. However, a generic feature of such 
models with a naked cusp in an inclined spiral is that the 
overall magnification is significantly lower than in large 
core lenses. Such a peculiar lens configuration, however, 
is expected to be rare. 

Lewis et al. (2002a) have developed a lensing model of 
the APM 08279+5255 system based on such a highly el- 
liptical lens, where the quasar core of the lensed galaxy is 
situated in the direct vincinity of a naked cusp. In their 
model, the inclined disk of the lens has a projected axis 
ratio of 0.25, and a rotational velocity of 200 kms^^. It 
has a core radius of 0.065 kpc, and is truncated at an 
outer radius of 8/i^^kpc. This truncated, flattened disk 
is assumed to sit in a spherical halo. Due to the possible 
detec tion o f microlensing (Lewis et al. 2002b; Chartas 
et al. |2002[ ). the observed optical and X-ray brightnesses 
of the individual images may not reflect the true rela- 
tive macrolensing magniflcations. Their brightness ra- 
tios thus are not considered a hard constraint for the 
lens mass model. Apart from the optical image posi- 
tions, Lewis et al. (|2002al ) also considered the structure 
seen in their C0( J=l— >0) map a constraint for the lens- 
ing model. While the continuum emission underlying the 
C0( J=l— i'O) emission in their maps is lensed into an arc 
and aligned with the optical emission, the CO(J=1^0) 
emission appears signiflcantly more extended. They thus 
assumed that the CO disk in APM 08279+5255 has to 
be large enough to reach out beyond the naked cusp 
and into the 5-image central caustic structure. The re- 
sulting lensing magnification of the optical emission is 
^opt_r,^ and image brightness ratios of 6(A,B)=0.75 and 
6(A,C)=0.75. The model reproduces the brightness ratio 
between images A and B (which is not subject to extreme 
changes) relatively well, while image C is significantly 
brighter than observed. The lensing magnification of the 
CO emission in their model is only /i2'-'=2.5-3. The true 
size of the CO disk thus is 400-1000 pc, and the CO be- 
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Fig. 12. — New lensing model for APM 08279+5255, assuming a highly inclined spiral lens. The top left panel shows the caustic and 
critical line structure. The orange circle indicates the source position, and the red circles represent the corresponding image locations. The 
right panel shows the image configuration for a range of source sizes. The color contours depict an increase in radius in 0.25 kpc steps. The 
bottom left panel shows the corresponding overall magnification as a function of radius. In this model, the CO(J=l— »0) emitting region 
has a true radius of ~0.5kpc or less. 



comes the dominant contributor to the total mass within 
the central kpc of the QSO. Due to the significantly lower 
lensing magnification factor compared to previous stud- 
ies, this model predicts that APM 08279+5255 is intrin- 
sically an extremely bright source at X-ray-to-FIR wave- 
lengths, rather than being modestly bright. 

Our new, improved high-resolution C0(J=1— >0) maps 
have r evealed that the extended CO arc found by Lewis 
et al. (|2002aD is a noise artifact, and that the lensed CO 
emission is coincident with the optical images. In con- 
trast to the underlying, slightly more extended contin- 
uum emission, it even shows the image structure seen in 
the optical. Assuming that the CO emissi on is more com- 
pact than assumed by Lewis et al. (|2002af ) and lies within 
the region of the naked cusp, their model describes these 
new observations quite well. Their model then suggests 
that the CO-emitting region has a true radius of 350 pc 
or less, and has a larger magnification factor close to that 
of the optical emission. However, this model was derived 
taking the extended CO structure found by Lewis et al. 
(|2002af ) as a constraint. To find the best-fitting model 
to existing observations, it thus is desirable to perform a 
new parameter study. 

5.3. Lens Modeling: New Model 

In light of our new observations of APM 08279-1-5255, 
we have conducted an extensive, systematic parameter 
study of the lensing configurations allowed within the 
constraints of its observed properties. All considered 
models assume lensing by an inclined spiral galaxy. 

As described above, multiwavelength observations of 
APM 08279-1-5255 at different epochs appear to indicate 
that the image ratios change, possibly due to microlens- 
ing. Thus, they cannot be considered hard constraints for 



the primary lens. Leaving the relative image brightnesses 
as free parameters, however, would sacrifice two obser- 
vational constraints, and, in fact, would under-constrain 
most lens models. Models considered good thus are re- 
quired to reproduce 6(A,B) relatively well (due to the 
relatively minor observed fluctuations), and to reproduce 
C as the faintest image. 

In Fig. [T21 we show the model that gives the over- 
all best fit to the observed properties at different wave- 
lengths within our study. The highly inclined disk of the 
lens in this model has a projected axis ratio of only 0.09. 
It is truncated at an outer radius of 15 h^^ kpc, and has 
a rotational velocity of 140 kms^^. The top left panel 
of Fig. [T2] shows the caustic and critical line structure, 
as well as the position of APM 08279-1-5255 in the source 
and image planes. The model predicts 6(A,B)=0.84, and 
6(A,C)=0.57, respectively. The brightness ratio of im- 
ages A and B is close to the observed value. The total 
magnification of a point-like source in this model is only 
4.15, and thus relatively modest. The bottom left panel 
shows the differential total magnification as a function of 
the source radius. The total magnification scales down 
by only a factor of 2 from source radii of few pc out to 
20 kpc, showing that differential lensing effects do not 
strongly influence the observed flux ratios. Out to 5 kpc, 
it varies by less than 20%. The SED shown in Fig. [H] 
thus is predicted to not be strongly distorted by differ- 
ential lensing effects over the whole spectral range. The 
right panel of Fig. [1^ shows the observed image con- 
figuration and overall morphology for a range of source 
sizes. The morphology of the C0(J=1— *0) observations 
is consistent with having a true source radius of 500 pc 
or less in this model. This is consistent with a lensing 
factor of 4 in Fig. [TI] assuming an inclination of <30° 
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and a high CO area filling factor. In the following, we 
will assume rtruc=550pc and ^l=4 for the CO(J=1^0)- 
emitting region, consistent with both observations and 
the lensing model. The model also correctly reproduces 
the 2.6 mm continuum emission underlying the CO line 
emission, assuming that the extended continuum emis- 
sion component comes from a region of about 1 kpc size. 
A more detailed description and analysis of the uncer- 
tainties within this model of the gravitational lensing in 
APM 08279+5255 will be described by G. Lewis et al. 
(2007, in prep.). 

6. BLACK HOLE, GAS AND DYNAMICAL MASSES 

6.1. Black Hole Mass 

The enormous energy output that powers the huge 
bolometric luminosity Lboi of APM 08279+5255 is pro- 
duced by a central super-massive black hole (SMBH) that 
is fed through an accretion disk. There are several ways 
to constrain the SMBH mass from observations. The 
most basic approach is to assume that the gradient of 
the radiation pressure driving gas out of the central re- 
gion is exactly counterbalanced by the gravitational force 
attracting infall of the same material. This equality of 
forces provides a limiting mass for the SMBH. Assuming 
that dominant mechanism contributing to the opa city is 
Thomson scattering off free electrons (Eddington Il921h 
leads to: 
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(9) 



Assuming Lboi^^Edd gives M||^'^= 2.1 x 10" ^^^M©. 

Another, perhaps better constrained estimate is pro- 
vided by observations of the 23^5*1/2-2^^^1/2,3/2 res- 
onance transitions of the C'^+ ion (C ivAA1548,1550). 
From reverberation mapping of nearb y active galactic 
nuclei, Vestergaard & Peterson ()2006D have found an 
empirical relation between the C ivAA1548,1550 FWHM 
linewidth Atip^^j^^, the rest-frame UV continuum lumi- 
nosity at 1350 A, and the mass of the central black hole: 



A/, 



BH 



Mf; 



4.6 X ( fi!taM 
kms 



ALA(1350A)y 
1044 erg s-i J 



(10) 

From the spectrum by Irwin et al. ()1998l ). we thus 
derive M§j|V=9.0 x 10^° hI^Mq, or about half of the 
Eddington limit. However, note that the reverberation- 
based masses underlying this equation are uncer tain by 
a factor of about 3 (Vestergaard & Peterson l2006l ). Also, 
part of the wings of the C iv line in APM 08279+5255 
are subject to strong absorption, making the derivation 
of the linewidth somewhat uncertain. However, the es- 
timated width of the C iv line is similar to that o f the 
Paa and Pa/3 lines of hydrogen (Soifer et al. I2004D . Fi- 
nally, the accretion disk may be inclined toward the line 
of sight, which would also influence this result. 

6.2. Molecular Gas Mass 

To derive molecular gas masses Afgas from CO 
luminosities in high rcdshift quasars, it currently 
is the most common approach to adopt a well- 
established conversion factor aco for nearby ULIRGs 



(Q^cq =0-8M(D(K kms~^pc^)~4; Downes & Solomon 
Il998l ). This is justified by the fact that the observed 
brightness temperatures in distant quasars appear to 
be similar to those of local ULIRGs. However, aco 
also depends on the gas density, which appears to be 
much higher in APM 08279+5255. Weifi et al. (j2007l ) 
thus re-derive a much higher conversion factor of aco 
~ 5 M0(K kms~^ pc^)~^ for this source. From our 
CO(J=1^0) luminosity, we thus derive Mgas=5.3 x 
10^^ Mq. This corresponds to about six times A/bh 
(based on the Civ estimate), assuming that no signifi- 
cant differential lensing effects act between the scales of 
the SMBH and of the molecular gas reservoir. This is 
somewhat low er tha n the factor of roughly ten found by 
Shields et al. (j2006f ) for a larger sample of distant gas- 
rich quasars, but in agreement within the large spread of 
observed values. 

6.3. Dynamical Mass 

From Newtonian mechanics, a dynamical mass 
can be derived for the rotating molecular disk in 
APM 08279+5255 (e.g., Solomon & Vandcn Bout 
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where Awp^jjjyj is the observed CO linewidth, and 
Ttiuo is the derived true radius of the CO disk. Using 
the observed C O(J=1 ^0) hnewidth of 556 ± 55 kms'^ 
(Riechers et al. I2006D . and using rtrue=550pc, we de- 
rive A^fdynSin^i = 4.0 x 10^" Mq. This dynamical mass 
of the central region of APM 08279+5255 corresponds 
to the sum of black hole mass, gas mass, dust mass 
and stellar mass (plus dark matte r) . A ssuming a gas- 
to-dust ratio of 150 (Weifi et al. 120071 ). the contribu- 
tion of dust to the total mass budget is negligible. For 
the moment, we will also assume that the mass frac- 
tion of the stellar bulge in the central 550 pc is small, 
i.e., that ATdyn sin^i~(Mgas+AfBH)- Assuming a face- 
on, unlensed disk, Afgas+AfBH would exceed the dynam- 
ical mass by a factor of 15.5, and thus is ruled out by 
the observations. Assuming /iL=4, we thus find that 
the CO disk in APM 08279+5255 can be inchncd by at 
most i=30° to not exceed the dynamical mass. With 
— A-,, CO sin~-'^z, this implies a rotational velocity 



AWrot^Aw- 

of 1095 kms 



FWHM 

-1 



6.4. Stellar Bulge Mass Limit 

As mentioned in the previous section, there is a 
third important constituent to the dynamical mass, 
i.e., the mass of the stellar bulge. This is a quan- 
tity that currently cannot be observed in the high- 
est redshift galaxies (particularly in quasars, where 
the AGN outshines any stellar light from the host). 
However, from our observations and theoretical con- 
siderations, we can derive a limit for the mass frac- 
tion o f the stellar bulge within rtruc- Carilli & Wang 
()2006l ) find that high-redshift gas-rich QSOs have a me- 
dian A'yp^jj,^ of 300 kms~^, at a mean inclination of 
(z)=13° (see their erratum). This corresponds to a 
mean rotational velocity of (Ai;i-ot)=1335 kms^^. As- 
suming that this value applies to APM 08279+5255, 
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TABLE 6 

Mass Budget in the Central Region of APM 08279+5255. 



Component 


Mass 




Mass Fraction 




lensed 


unlensed 






[IQio (/.Z')M0] 


[IOIOMq] 


[%] 


Molecular Gas 


53 


13 


57.4 


Dust 


0.35 


0.09 


0.4 


Black Hole 


9.0 


2.3 


9.7 


Stellar Bulge'' 


30 


7.5 


32.5 



Note. — Here, wc assume i— 25°, rtruc— 550pc, and iil—A. 
Masses are uneorrected/corrected for lensing, assuming 
^ Fraction within the radius of the CO disk. 



we derive an inclination of i=25°. Assuming 
MdynSin2i=(Mgas+Afdust+A^BH+A4uigo), and a gas-to- 
dust ratio of 150, this corresponds to a bulge mass 
of Mbuige=3.0 X 10"m7^Mq, or Mbuige=3.4MBH 
(0.57 Mgas; see also Tab. |6| for /xl=4. For comparison, 
if a model with /zl = 100 (and thus rtrue=125pc) that 
fits the data were to exist, it would give Mbuige=4.4 x 

10^2 Mq, or Mbuige=50MBH (8.4 Mgas). This means 
that for low-/iL models, the gas component dominates 
the total baryonic mass budget in the central region of 
APM 08279-^5255, while for high-/XL models, the stellar 
bulge is the dominant contributor. 

6.5. Limitations of this Analysis 

The maximum CO velocity widths predicted by the 
Carilli & Wang (|2006[ ) unified model exceed the rota- 
tional velocities of the most rapidly rotating nearby giant 
ellipticals by about a factor of 2. This may suggest that 
the observed linewidths in the systems with broader but 
not multiply peaked CO lines are not due to rotationally 
supported disks, but possibly in a later stage of a major 
merger where the disk is not self-gravitating yet. In that 
case, the mean inclination of quasar host galaxies derived 
by Carilli & Wang may overpredict the rotational veloci- 
ties for gravitationally bound molecular disks. This is of 
particular concern for APM 08279-1-5255, as its CO ve- 
locity width is already the highest observed among the 
distant quasars. If the measured CO linewidth was that 
of a gravitationally bound disk and representative for 
the rotational velocity inside a present day giant ellipti- 
cal, the inclination could not be much lower than z=60°. 
This leaves a factor of 2-3 discrepancy in the Mdyn esti- 
mate relative to the i=30° solution. 

The uncertainties in the luminosity-based mass esti- 
mates are a factor of a few. Together with the uncertain- 
ties in disk size and lensing factor, and the contribution 
of dark matter to the mass budget, the real inclination 
of the molecular disk remains uncertain (as only a factor 
of 3-4 difference exists between high and low inclination 
solutions). Due to the degeneracy of area filling factor 
Ja and inclination, our r esults thus are still in agreement 
with those of Weifi et al. (|2007f ) within conservative error 
estimates, even with the new morphological constraints. 

7. ON THE RADIO-FAR-INFRARED CORRELATION 

In the local universe, star-forming galaxies follow a 
tight correlation between the monochromatic radio con- 
tinuum luminosity at 1.4 GHz and the far-infrared lumi- 
nosity (ipiR): which holds ov er mo re than 4 orders of 
magnitude (e.g. , Con don et al. Il99lt Condon 1992; Yun 
et al. 120011: Bell boOSD . This correlation is probably due 



to coupling between infrared thermal dust emission due 
to heating by young stars and non-thermal radio syn- 
chrotron emission from supernova explosions and SNRs 
(e.g., Helou etal. |1985f ), and constitutes the basis for 
the Condon model described above. The radio-FIR cor- 
relation can be ex pressed via the so-called q-parameter 
(Helou etal. Hill): 

'^ = ^°H 9-8xlO-i5Lj '^"HwHF^j- (''^ 

The apparent FIR luminosity of APM 082 79-^5 255 is 
Lftr =(2.0 ±0.5) X ^0^'^^il^LQ (Beelenetal.UMi Weifi 
et al. l2007l ). We thus obtain g=2.03, a value which falls 
well within the range of values found for IRAS galaxies 
(Yun etal. 1200 ll ). and consistent with the values found 
for o ther dust and gas-r ich high-z q uasars (CarilH etal. 
120041 : Beelen et al. l2006t Wang et al. l2007| ). A common 
interpretation of a g value close to the radio-FIR corre- 
lation is that star formation is the dominant process to 
power both the radio and FIR emission. However, us- 
ing the dust model of Weifi etal. (j2007f ). our analysis of 
APM 08279-^5255 has shown that the AGN may power 
the dominant fraction of both the FIR and the radio 
emission. We obtain this conclusion based on the fact 
that the SED of APM 08279+5255 cannot be fitted with 
a simple Condon model (which motivates the star for- 
mation interpretation of the radio-FIR correlation). We 
thus conclude that caution has to be used when inter- 
preting the properties of dust and gas-rich high-z quasars 
based on the radio-FIR correlation alone. 

8. DISCUSSION 

The observed properties of APM 08279-^5255 are un- 
doubtedly peculiar. It thus is important to understand 
whether this peculiarity is due to a rare lens configu- 
ration, possibly leading to very high lensing magnifica- 
tion factors and significant differential lensing effects, or 
whether APM 08279-1-5255 is a source that is just mod- 
erately magnified by gravitational lensing and thus in- 
trinsically extreme. 

8.1. Multiwavelength Properties of APM 08279+5255 

APM 08279-1-5255 shows broad absorption lines in the 
optical/UV, and emission from high excitation lines like 
Civ and Nv close to the central nucleus (Irwin etal. 
Il998f ). kinematically blueshifted relative to the molecu- 
lar g as and Ci emissio n by about 2500 kms~^ (Downes 
et al. fl999t Wagg et al. [2001 . The quasar also shows rel- 
ativistic X-ray broad absorption lines from even higher 
emission lines of iron from gas situated within the UV 
BAL region (Chartas etal. 2002; Hasinger etal. 1200^ . 
The BALs come from an outflow of highly ionized gas 
from the accretion disk, mostly driven by radiation pres- 
sure. This outflow distributes the accretion disk material 
over the central region of the quasar, and out into the 
host galaxy. The X-ray BAL region is likely a dust-free, 
high column density absorber responsible for shielding, 
as indicated by the detecti on of an iron K-shell absorp- 
tion edge (Hasinger etal. l2002f ). This shielding is re- 
sponsible for the typical X-ray faintness of BAL QSOs. 
The X-ray luminosity alone thus is not a good measure 
for the energy output of APM 08279-f 5255. The detec- 
tion of strong iron lines also indicates that the X-ray 
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BAL emission comes from significantly metal-enriched 
gas with possibly super-solar metallicities. While this 
strongly indicates that the gas has already been pro- 
cessed in a starburst cycle, it is biased toward the very 
central region (possibly <0.1pc). It thus cannot be as- 
sumed ad hoc that the metallicity is super-solar over the 
whole scale of the galaxy. The continuum X-ray emis- 
sion in APM 08279-1-5255 has a photon index that is in 
qualitative agreement with a radio-quiet quasar (RQQ). 
The optical/IR shape of the SED is also consistent with 
APM08279-H5255 being a RQQ. Even more so, it fol- 
lows the radio-quiet locus of the radio-FIR correlation, 
which indicates that the copious amounts of dust in this 
source may be heated to a significant fraction by newly 
formed young stars, while the radio synchrotron emis- 
sion may originate from supernova remnants (Irwin et al. 
Il998t Beelen et al. |2006[) . However, our analysis of the 
radio and dust properties indicates that most of the ra- 
dio/FIR emission may well be AGN- related. Within the 
unified theory of AGN galaxies, the smoothness of the 
IR spectrum, in particular the lack of a pronounced IR 
peak, indicates that the accretion disk in thi s radi o-quiet 
quasar is seen close to face-on (Soifer etal. l2004[ ). This 
agrees well with the enormous optical brigh tness and the 
flatness of the optical spectrum (Irwin et al. '1998'; Egami 
etal. HOOO), indicating a lack of obscuration and thus 
an angle of view which clearly lies within the ionization 
cone. The hot dust in the parsec-scale accretion disk is 
heated by the central AGN, and extends outwards into 
a warm dusty torus, extending out to 100 pc scales. 

The luminous AGN in APM 08279+5255 is the dom- 
inant power source out to large scales, which even may 
generate a large fraction of the huge (far-)infrared lumi- 
nosity in this source by heating the dust.^* The appar- 
ent FIR luminosity of APM 08279+5255 is Lfir=(2^± 
0.5) X 10^^ /i7^L w (Beelen etal. ^2006; Weifi etaL [200l . 
Weifi et al. (|2007l ) have found that about 90% of the FIR 
luminosity are due to emission from warm 220 K dust 
which may be powered by the AGN, while only 10% 
are due to colder 65 K dust that is likely heated by star 
formation. This would explain why a significant frac- 
tion of the dust in APM 08279+5255 on extended (100s 
pc) scales is warmer than in typical ULIRGs, or even 
other high-z, dust-rich QSOs. It would also explain why 
APM 08279+5255 is an outher on the L'^^^q-Lfiyi relation 
(Riechers et al. 2006) , which is set by galaxies where the 
FIR luminosity is dominated by heating from star forma- 
tion. Subtracting the 90% of Lfir that may be powered 
by the AGN, APM 08279+5255 follows the L'^q-Lfir re- 
lation quite well. However, note that even though the ap- 
parent FIR luminosity of APM 08279+5255 is extreme, 
this would indicate that only 3% of the total quasar 
power (as measured by Lboi) are re-emitted in the FIR 
wavelength regime (assuming no differential lensing ef- 
fects play into the ratio). We find that up to 80-90% of 
the radio luminosity in APM 08279+5255 may be pow- 
ered by the AGN, and only 10-20% by the starburst. 
Due to the fact that this ratio is similar to that found 
for the FIR luminosity, APM 08279+5255 stiU follows the 
radio-FIR correlation for star-forming galaxies. 

^® It is important to note that, due to the unique properties of 
APM 08279+5255, additional diagnostics are desirable to undoubt- 
edly disentangle the heating sources of the gas and dust. 



APM 08279+5255 thus is a radio-quiet, optically lu- 
minous dust-rich quasar in which the central accretion 
disk and dust torus are likely seen close to face-on (note 
however the above concern about the inclination). Its 
high FIR luminosity can be explained by dust heating of 
the AGN, as can be the high temperature of the dust. It 
even explains the high CO excitation if the molecular gas 
is situated in a rather compact (and thus dense) circum- 
nuclear ring seen at a low inclination toward the line of 
sight, just like the central region. The high dust and gas 
temperatures and relatively high densities also explain 
the exceptional HCN (and HCO"*") excitation in this sys- 
tem, as IR-pumping via the V2 — I vibrational bending 
mode beco mes efficient at such high temperatures (Weifi 
etal. I2OO7D . Without shielding and self-shielding, the 
gas and dust would be even warmer than observed. As- 
suming that only the fraction of Lboi re-radiated in the 
FIR reaches out to the molecular regions, the temper- 
ature at a radius of 125 pc would still be 200 K. Even 
at 550 pc, the temperature would not drop significantly 
below 100 K. The LVG-predicted temperatures for the 
cool, dense gas component (Tkin — 65 K) thus require 
self-screening if the gas is indeed located in a 550 pc re- 
gion around the AGN.^^ Note that all the above may be 
considered extreme (i.e., rarely observed). However, all 
of these effects can be explained physically without as- 
suming extreme and/or differential lensing effects. The 
CO luminosity in APM 08279+5255 is relatively modest 
(compared to other high-z sources) even without correct- 
ing for lensing. If Lfir is corrected assuming a lensing 
magnification factor of only 4, as suggested by our new 
model, it is not the highest observed value even with- 
out correcting for AGN heating. Both Lqq and LpiR 
are higher for the z=4.7 quasar BR 1202-0725, which is 
thought to be unlensed, but under going a massive, gas- 
rich merger event^° (Carilli et al. l2002t Riechers et al. 
[2006h . 

8.2. The Black Hole in APM 08279+5255 

Correcting for /xl=4 predicts A/bh=2.3 x 10^" M© for 
APM 08279+5255. This appears to be extreme com- 
pared to more typical SMBH masses of a few times 
10^ Mq of high-z quasars. However, comparable or even 
higher black hole masses are found for BRI 1335-0417 
(z=4.41) and BR 1202-0725 (however, the estimates for 
these sources are based on strongly absorbed C iv emis- 
sion lines; Storrie-Lombardi etal. 1996i : Shields etal. 
I2OO6I ). Moreover, Vestergaard (20(3has found that 
for the range of black hole masses, an upper enve- 
lope of Mbr^IO^^Mq is observed over a large range 
in redshift (4<z<6 for her sample). So, again, while 
APM 08279+5255 appears relatively extreme, even its 
AGN properties are far from being unique in a sce- 
nario assuming only /Xi=4. The extreme optical luminos- 
ity may then be explained by a luminous AGN episode 
of such a few times 10^° Mq SMBH, which accretes at 
(super-)Eddington rates (within the uncertainties of the 
BH mass estimators). Also, the formation of such a 

Due to the fact that T oc /i^"'^''* and ''true Ml/=4 

model requires more screening than the previous /ii=100 models. 

20 Both the CO and FIR emission in BR 1202-0725 are however 
distributed over at least two distinct components, which then each 
have comparable but lower luminosities than APM 08279+5255. 
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SMBH at the high mass end of the BH mass function 
by a redshift of 4 is consistent with 'downsizing', which 
predicts that such high mass black holes form at the high 
density peaks at high redshift, and are built jup^ rapidly 
by vigorous accretion (e.g., Di Matteo etal. 2007). Due 
to the fact that the SMBH in APM 08279+5255 is among 
the most massive black holes that are observed, it is likely 
to end up as a central cD galaxy of a massive cluster. 
For such galaxies, recent cosmological simulations sug- 
gest that a large fraction of the mass of the black hole 
grows from mergers with other black holes rather than 
bein g accr eted by a single progenitor black hole (Sijacki 
etal. IMS). 

8.3. The Stellar Bulge of APM 08279+5255 

In addition to Mbh, we were also able to derive 
a limit on the mass of the stellar bulge (Mbuigo) for 
APM 08279+5255. However, our dynamical mass esti- 
mate, which is used as the predictor for Mbuige, is only 
valid for the inner region of the galaxy where the bulk of 
the molecular gas is found. For the /iL=4 model, this de- 
scribes a region with a radius of 550 pc. For galaxies with 
^^BH>10^ Mq used in local studies of the MeH-crbuigc re- 
lation, it has been found that the host galaxies are typ- 
ically giant elliptical galaxies with effe ctive radii of 1.5- 
8kpc (Tremaine etal. 2002; Faber et al. ll997| ). Assuming 
that the volume density of stellar luminosity (and thus 
mass) flattens with radius from to in th e inne r 
region of the elliptical host galaxy (Gebhardt et al. ll996f ). 
or even that it flattens into a 'core', a significant fraction 
of the bulge mass is expected to be found inside the in- 
ner 550 pc. From the local AfBH~-^^buigc relationship, it is 
expected that Mbuigc'~700MBH- For the inner region of 
APM 08279+5255, we have found Mbuige=3.4MBH- Al- 
though uncertain within a factor of a few, this value falls 
short by more than two orders of magnitude compared to 
the local estimate. Even if corrected for a lOkpc radius 
host galaxy, a clear offset from the local MBH--Mbuigc 
relationship remains^ ^. This is consistent with what is 
found for other quasars at even higher redshift studied 
with_thc s ame tec hnique (Walter et al. [2004; Riechers et 
al. l2008ai l2008bl : see also discussions by Shields etal. 
120061: Ho [2007fr Based on different techniques, similar 
trends are found for various galaxy s amples at low and in- 
termediate reds hift (e .g., Peng et al. 12006 1; McLure et al. 
120061 Woo et al. l2006t Salviander et al. :2007:). This may 
indicate a breakdown of the local MBH-Afbuigc relation- 
ship toward high redshift, and that the mass correlation 
between SMBHs and their stellar bulges is not a uni- 
versal property, but rather an endpoint of a long-lasting 
evolutionary process throughout cosmic times. We thus 
conclude that the SMBH in APM 08279+5255 appears 
to already be largely in place, while the buildup of the 
stellar bulge is still in progress. 

8.4. A Possible Evolutionary Scenario for 
APM 08279+ 5255 

If the assumptions made to derive the properties of 
APM 08279+5255 are correct, this result has several in- 
teresting implications. It appears as if SMBH growth 

Note that the local AfBH-fbulgc relationship is difficult to 
reconcile even with a fj,j^=100 model, which leaves more room for 
Mbuigc within M^yn (i.e., Mbuige=50 Mbh)i but also has to assume 
host galaxy size relations for less massive black holes. 



and formation through gas dissipation and accretion (and 
possibly mergers) was the main characteristic during the 
early formation of this source (even relative to star for- 
mation). Since the SMBH in APM 08279+5255 is al- 
ready found at the high mass end of the black hole mass 
function, it is expected that most of the molecular gas 
present at z=3.9 is either blown out or formed into stars 
by z—0, while only a minor fraction is accreted onto the 
black hole. Moreover, if APM 08279+5255 is to evolve 
into a galaxy that fulfills the AfBH"0'buigc relation by 
z—0, a significant fraction of the molecular gas has to ac- 
tually end up in stars. In fact, even if the whole amount 
of molecular gas that is currently present in the host 
galaxy of APM 08279+5255 were to be converted into 
stars by z=0, it would not be sufficient to reach the lo- 
cal MBH-cTbuigo relation. If our picture of this galaxy is 
correct, this result implies that the buildup of the stel- 
lar bulge cannot be accomplished by only converting the 
observed molecular gas reservoir into stars, but relies to 
a significant part on the accumulation of stellar matter 
via other mechanisms, and/or the the accretion of ad- 
ditional star-forming material. Although further details 
are difficult to quantify at present, it thus is likely that 
a significant fraction of the spheroid will be produced 
via major and minor mergers (which are also required 
for the re-distribution of angular momentum). It has 
been suggested that a common, growth-limiting feedback 
mechanism acting on both star formation and black hole 
assembly is the physical explanation for the local Mbk- 
fbuigc relation. Assuming that such a mechanism exists, 
quasar winds as a form of AGN feedback would be an ob- 
vious c andid ate for a system like APM 08279+5255 (Silk 
& ReeslHH). One difficulty with this assumption would 
then be that the SMBH in this source already accretes 
radiatively highly efficiently at (super-)Eddington rates 
(see above), where such feedback effects to regulate the 
further growth of the SMBH are expected to be very 
strong. If this effect were to also regulate star formation 
in the host galaxy, it would be expected to shut off both 
on a relatively short timescale by pushing the remain- 
ing gas outwards. The build-up of the stellar bulge then 
can only proceed further through hierarchical merging 
and/or accretion of additional gas, which is expected to 
again feed both the SMBH and star formation. To end 
up anywhere near the local MBH~crbuigc relation through 
a self-regulated mechanism, APM 08279+5255 thus is re- 
quired to form and/or accumulate stars more efficiently 
(relative to the black hole feeding/merger rate) than in 
the past. One possibility to build up stellar mass while 
avoiding significant SMBH growth ma y be t hrough sub- 
sequent 'dry' mergers (e.g., Rines et al. l2007t ). which may 
act in combination with other processes to assemble the 
massive host of APM 08279+5255 by z=0. 

8.5. Our Current Picture of APM 08279+5255 

There are two main points that support the iil=^ lens- 
ing model proposed for APM 08279+5255. First, the 
observed CO(J=1^0) properties give rather tight con- 
straints on the allowed combinations of intrinsic CO disk 
size and magnification strength for the favoured range 
of low disk inclinations, assuming a relatively high CO 
area filling factor. The /iL=4 model lies well within the 
range of these constraints. Second, the spatially resolved 
multiwavelength observations of APM 08279+5255 show 
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similar morphologies of the lensed source from the X-ray 
to radio wavelength regime, probing sub-pc to hundreds 
of pc scales, while showing a different morphology out 
to kpc scales. The /iL=4 model produces this quite nat- 
urally, predicting modest effects of differential lensing. 
The high lensing magnification models suggested so far 
do not reproduce the second characteristic. Also, for 
low inclinations, they predict a very compact CO disk of 
only about 100 pc radius. This is by a factor of a fe w less 
than in typical ULIRGs (Downes & Solomon 1998). The 
/iL=4 model, on the other hand, produces a CO disk of 
550 pc radius, which is in a more typical regime. 

While some properties of APM 08279-^5255, like the 
optical luminosity and expected strong tidal forces from 
the black hole, as well as the requirement of significant 
shielding of the molecular gas, may still favour a sce- 
nario in which the source is highly magnified by grav- 
itational lensing, there thus are several arguments that 
favour models with low /i^- It however is important to 
note that, if a high-/ii model was to be found that fits the 
observed morp hologi cal properties of APM 08279-f 5255 
(see Krips et al. l2007l for a recent attempt), further stud- 
ies would be necessary to exclude one scenario or the 
other. Based on the existing observations, no definitive 
conclusion can be drawn - the final decision can only 
be made by detecting the lens itself. The lensing model 
presented in this paper aids in developing a quite consis- 
tent picture of APM 08279-1-5255 based on the existing 
observations, describing the source as a dust and gas- 
rich galaxy with a very massive, active black hole. The 
deep gravitational potential of the central region of this 
source causes the molecular gas to be rather dense, while 
the strong, penetrating AGN radiation causes it (and 
the dust) to be rather warm. The whole system may be 
seen relatively close to face-on (the inclination is however 
difficult to constrain). It has a co-eval starburst which 
only contributes about 10% to the FIR dust heatin g, but 
still produces more than 500MQyr-i (WeiB etal. l2007l . 
correcting for ^il=A). The star formation in this galaxy 
takes place far off the black hole mass-bulge mass scal- 
ing relation of nearby galaxies, suggesting that the latter 
may be subject to significant evolution throughout cos- 
mic times. 

9. SUMMARY 

In this paper, we present improved, high resolution 
C0( J=l— *0) observations toward the z=3.9 BAL quasar 
APM 08279-f 5255 with the VLA. We also present multi- 
wavelength radio continuum observations, and a search 
for CN(iV=1^0) emission. The source's properties are 
investigated based on a revised model of the gravitational 
lensing in this system. 

From our comprehensive analysis, we obtain the fol- 
lowing main results: 

1. The CO(J=1^0) emission is resolved and has a 
(lensed) source size of ~0.3". The structure shows a mor- 
phology that is very similar to the X-ray/optical/NIR 
morphology, indicative of three images of a compact, 
gravitationally lensed source (of which two are blended 
in the CO maps). In particular, we do not find any ev- 
idence for extended CO(J=1^0) emission over a scale 
of several arcseconds, iiicontrast to previous reports 
(cf. Papadopoulos et al. [200lh . The whol e sing le-dish 
CO(J=1^0) flux found by Riechers et al. (|2006f ) is re- 



covered from the compact structure detected in our VLA 
maps. The rest-frame 2.6 mm continuum emission un- 
derlying the CO(J=1^0) line is slightly more extended 
than the CO source, and may correspond to differen- 
tially lensed emission either due to free-free radiation 
from the starburst or extended AGN-related emission (in 
agreement with the overall steep spectral index). A ra- 
dio source is detected 3" south of APM 08279-^5255 in 
two out o f six radio continuum bands (see also Ibata 
et al. ll99"9l ). but probably not related to the high redshift 
quasar. 

2. We derive an apparent CO(J=1^0) luminosity of 
^co(i-o) = (10-6±0.9) X lO^VZ^ K kms-ipc2. This is 
by about 20% lower than the re-derived C0( J=2— >1) lu- 
minosity, and even by 40% lower than the CO(J=4^3) 
luminosity. This may indicate that the optical depth 
in the lower J CO transitions is only modest, as al- 
ready indic ated b y previous studies (Riechers et al. 120061 : 
Weifi etal. l2007f) . Due to the high gas density, the 
C0( J=1^0) luminosity corresponds to a total apparent 
II2 molecular gas mass of Mgas=5.3 x 10^^ /x^^ M0. 

3. We set a limit on the CN(A^=1^0) luminosity in 
this source (based on GBT spectroscopy), which is <37% 
of the CO(J=1->0) luminosity. 

4. The radio continuum SED of APM 08279-^5255 is 
likely dominated by AG N emis sion. Based on a starburst 
galaxy model (Condon Il992f ). we conclude that free- 
free and synchrotron emission from the host galaxy con- 
tribute <20% to the radio continuum emission (note that 
this fraction is clearly model-dependent). We also find 
that the AGN lies in the transition region between radio- 
quiet and radio-loud quasars, and is in overall agreement 
with being radio-quiet. The low frequency part of the 
radio SED is (ultra-)steep, which is commonly found in 
high-z quasars with extended radio emission. The over- 
all spectral slope of the radio SED is in agreement with 
synchrotron emission that may originate from both the 
starburst and the AGN. 

5. In spite of the fact that both the FIR (WeiB etal. 
I2007t) and radio continuum emission in APM 08279-1-5255 
appear to be dominantly powered by the AGN, it follows 
the radio-FIR correlation for star-forming galaxies with- 
out correcting for the AGN contribution to the radio/FIR 
luminosities (^80-90% in both cases). This may indicate 
that loyalty to the radio-FIR correlation does not nec- 
essarily imply that a galaxy's emission is star formation 
dominated on an object-to-object basis. 

6. From our CO maps, we derive observed brightness 
temperatures of y°'''^^^+^~23 K for the blend of optical 
quasar images A-fC and r^'^'''^~13K for quasar image 
B. Assuming a true brightness temperature of 50 K as 
predicted by LVG models of the CO excitation (Weifi 
et al. I2007f) . this corresponds to a lensed CO equivalent 
disk with a radius of ~900pc. 

7. Our revised lensing model of APM 08279-1-5255 pre- 
dicts that the luminosity is only enhanced by a lensing 
magnification factor of about /iL=4 over the whole ob- 
served wavelength range. Assuming that the CO disk is 
seen relatively close to face-on (i<30°) and has a large 
area filling factor, this model is consistent with a CO ring 
that has an unlensed radius of about 550 pc. 

8. From C iv observations of APM 08279-f 5255 and 
a scaling relation of nearby galaxies from the literature. 
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we derive an apparent black hole mass of A/bh=9.0 x 
10^° Mq, which corresponds to about half of the Ed- 
dington limit. This also corresponds to about 17% of the 
molecular gas mass. 

9. From our CO maps and model, we derive a 
dynamical mass of MdynSin^i = 4.0 x 10"^° Mq for 
APM 08279+5255. Assuming a<l=4, this indicates an 
inclination of i<30°, consistent with observations at all 
wavelengths within a (limited) unified scheme. Assum- 
ing an average rotational velocity for high- 2; quasars, this 
even sets a limit on the mass of the stellar bulge within 
the central 550 pc of Mbuige=3.4MBH, which falls by al- 
most two orders of magnitude short of the MBH--^buigo 
relation obse rved fo r nearby galaxies (see also discussion 
by WeiB et al. l2007f) . We thus conclude that the SMBH in 
this distant quasar is already in place, while the buildup 
of the stellar bulge is still in progress. A similar indica- 
tion was found for the z=6.4 quasar SDSS J1148-f 5251 
(Walter etal. 2004), the z=4.1 quasar PSS J2322+1944 
(Riechers et al. ,2008a'). and the z—AA quasar BRI1335- 
0417 (Riechers et al. 2008b), and thus suggests a break- 
down of this relation at early cosmic times and/or toward 
the high mass end. 

The observations presented herein highlight the im- 
portance of spatially and dynamically resolved studies 
of molecular gas in galaxies at z^A for our general pic- 
ture of galaxy formation and evolution. In the next 
decade, the Expanded Very Large Array (EVLA) will 
enable us to commence such studies at improved sen- 
sitivities and spectral resolution. The low J transition 
studies of molecular gas with the EVLA will be comple- 
mented by observations of the higher J transitions with 
the Atacama Large (sub-)Millimeter Array (ALMA) at 
even higher spatial resolution. Both instruments will al- 



low us to probe significantly deeper down the luminosity 
function(s) of the different populations of high redshift 
galaxies using molecular gas studies, and thus largely im- 
prove our picture of galaxies in the early universe. 
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